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INTRODUCTION 


Combustion  in  liquid  rocket  motors,  like  that  in  the  diesel  engines  or  gas  turbines, 
occur  as  a  rule  under  fully  turbulent  hydrodynamics.  The  turbulent  fluctuations  lead 
into  complex  flow/flame/droplet  interactions  which  substantially  modify  the  droplet 
evaporation  and  combustion. processes  (Faeth,  1983).  Because  of  the  large  number  of 
parameters  that  influence  the  reactive  flow,  systematic  investigation  of  turbulent 
spray  combustion  requires  consideration  of  small  local  flow  models  which  aim  at 
illuminating  one  particular  aspect  of  the  complex  problem.  For  example,  many  model 
studies  consider  combustion  of  a  single  droplet  or  an  array  of  droplets  in  stagnant 
atmosphere  or  under  specified  flow  conditions.  The  results  of  such  studies  could  then 
be  collectively  employed  in  the  development  of  more  comprehensive  models  of  the 
turbulent  spray  combustion. 

In  the  following,  a  general  description  of  the  results  of  the  research  studies  aimed 
at  the  understanding  of  the  various  aspects  of  turbulent  spray  combustion  will  be 
presented.  In  particular  ,  the  relevance  of  the  model  studies  to  the  simulation  of  the 
local  flow/flame  conditions  within  turbulent  combustion  in  liquid  rocket  motors  will  be 
emphasized.  The  detailed  discussion  of  each  study  are  presented  in  the  Appendices  1 
through  VI  at  the  end  of  the  present  report.  In  the  sequel,  a  brief  description  of  major 
findings  in  each  of  the  studies  will  be  presented,and  their  impact  on  the  global  aspects 
of  turbulent  spray  combustion  will  be  assessed.  In  addition,  some  suggestions  for 
research  areas  in  need  of  further  future  exploration  will  be  identified.  Also,  some 
discussions  will  be  presented  as  to  how  the  ideas  and  the  results  obtained  in  these 
investigations  may  help  the  modeling  of  turbulent  spray  combustion  in  liquid  rocket 
engines. 


SPRAY  COMBUSTION  IN  ST  AGNATION-POINT  FLOW 

The  detailed  discussion  of  the  study  on  the  combustion  of  sprays  in  the 
stagnation-point  flow  has  been  presented  in  Appendices  I  and  II.  As  is  described  in 
the  Appendix  I,  turbulent  spray  is  considered  as  a  jet  of  liquid  hydrocarbon  spray 
which  is  injected  into  a  turbulent  oxidizing  environment,  a  condition  which  closely 
simulates  rocket  motors  using  hydrocarbon  fuel.  Therefore,  from  the  central  core  of  the 
spray  to  the  far  field  regions  of  the  oxidizing  atmosphere,  an  entire  spectrum  of 
concentrations  from,  pure  fuel  (F),  fuel  -rich  (F+eO),  near-stoichiometric  (F+O),  fuel- 
lean  (eF+O),  and  pure  oxidizer  (O)  may  be  encountered.  Such  regions  of  gaseous 
mixtures  will  be  under  constant  random  motion  and  may  contain  fuel  droplets. 
Therefore,  the  concentrations  of  the  regions  will  be  changing  because  of  the  turbulent 
mixing,  molecular  diffusion  and  droplet  evaporation. 

In  the  presence  of  chemical  reactions,  the  interfaces  between  the  (F)  and  (O) 
regions  identified  above  may  support  diffusion  flames.  Also,  nartiallv  premixed  flames 
will  be  encountered  between  the  regions  (F)-(eF-t-O)  and  (0)-(F+  eO)  with  the  flame 
receiving  either  fuel  or  oxidizer  from  both  sides  of  its  reaction  zone.  Finally,  premixed 
flames  will  occjr  in  regions  which  have  near  stoichiometric  concentrations.  It  is 
emphasized  that  while  the  premixed  flames  are  dynamic  because  of  their  propagation. 


the  diffusion  flames  only  move  due  to  the  motion  of  the  burning  droplets  or  the  bulk 
convection  of  the  burning  (F)-iO)  interfaces.  Under  the  velocity  fluctuations,  the 
interactions  between  the  droplets  and  the  randomly  moving  flame  surfaces  are 
inevitable.  Clearly,  the  understanding  of  the  problem  requires  knowledge  of  the 
transient  interactions  between  the  flow,  the  flame,  and  the  droplets. 

One  of  die  important  problems  to  be  addressed  in  such  turbulent  spray 
combustion  is  the  behavior  of  the  individual  droplets  as  they  approach  and  pass 
through  the  flame  surfaces.  Because  of  the  expansion  of  the  gas,  the  velocity  normal 
to  the  flame  surface  will  experience  a  sudden  jump  across  the  flame.  Another 
significant  hydrodynamic  parameter  which  influences  the  flame  behavior  is  the  rate  of 
change  of  the  flame  surface  area,  i.e.  flame  stretch  (Williams,  1985).  It  is  expected 
that  the  rate  of  strain  will  also  have  substantial  effects  on  the  spray  combustion 
similar  to  that  in  homogeneous  turbulent  combustion.  We  emphasize  that  thus  far 
little  is  known  about  the  subject  of  droplet-flame  interactions,  which  was  one  of  the 
fact  that  motivated  the  present  research  effort.  In  addition,  the  influence  of 
aerodynamic  stretching  on  the  behavior  of  flames  within  sprays  is  almost  totally 
unknown. 

In  the  Appendix  I,  the  results  of  our  study  on  the  behavior  of  reactive  and  non¬ 
reactive  sprays  in  stagnation-point  flow  are  discussed.  An  important  finding  in  these 
tests  was  the  pronounced  change  in  the  droplet  slip  velocity  as  it  crosses  the  flame 
surface.  Thus,  the  expansion  of  the  gas  within  flames  resulted  in  appreciable  droplet 
deceleration  (acceleration)  in  the  pre-flame  (post-flame)  regions.  Clearly,  the  change 
in  the  droplet  slip  velocity  will  have  substantial  effects  on  the  droplet  evaporation  and 
combustion.  We  emphasize  that  the  knowledge  of  the  droplet  trajectory  as  it  passes 
through  flames  are  quite  significant  for  our  understanding  of  turbulent  spray 
combustion.  Possible  application  of  the  experimental  technique  for  the  evaluation  of 
flame  propagation  velocity  in  sprays  was  also  described  (Chen,  et  al.,  1988).  Finally, 
the  addition  of  water  spray  to  lean  methane-air  flames  was  shown  to  result  in  the 
formation  of  a  yellow-orange  radiation  layer  downstream  of  the  flame  surface.  We 
note  here  that  the  problem  of  soot  formation  and  accumulation  poses  severe  difficulties 
in  methane  and  propane- fueled  rocket  motors.  More  detailed  description  of  the 
findings  are  available  in  the  Appendix  I  (Chen,  et  al..  1988). 

Another  important  result  of  the  experimental  investigation  is  concerning  the 
occurrence  of  acoustic  instabilities  in  the  sprays  when  certain  critical  values  of  the  fuel 
concentration  or  the  flow  velocity  are  exceeded.  In  particular,  the  concentration 
domains  for  acoustic  flame  were  determined  for  both  methane-air  as  well  as  methane  - 
air-ethanol  spray  systems,  as  described  in  the  Appendix  II.  The  coupling  with  the 
acoustic  waves  produced  within  the  burner  cavity  were  found  responsible  for  the 
observed  phenomenon  (Zinn,  1986).  Because  the  present  experimental  model 
provides  steady  and  stationary  flames,  the  diagnostics  of  the  properties  and  structure 
of  such  acoustically  unstable  flames  is  facilitated.  Another  significant  finding  in  the 
experimental  investigation  was  the  occurrence  of  the  so  called  hysteresis  phenomenon 
shown  in  Fig.  1.  In  the  tests,  the  contour  of  maximum  fuel  molar  concentrations  Xp  at 
various  nozzle  velocities  V  corresponding  to  the  onset  of  acoustic  flame  instability 
were  determined.  Next,  the  fuel  concentration  of  the  already  acoustic  flames  was 
gradually  reduced  until  the  flames  ceased  to  be  acoustic  and  the  minimum  critical  fuel 
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concentration  was  noted.  As  shown  in  Fig.  1 ,  the  two  contours  for  the  critical  values  of 
fuel  concentrations  just  defined  do  not  coincide  leaving  a  hysteresis  band  in  the 
middle.  This  behavior  is  a  direct  indication  of  the  significance  of  the  role  of  the  non¬ 
linear  phenomenon  in  the  process  of  acoustic-chemically  coupled  flame  instabilities 
(Zinn,  1986  ).  Clearly,  the  results  have  direct  bearing  on  the  understanding  of  the 
mechanisms  foi  initiation/suppression  of  acoustic  instabilities  and  is  in  need  of  further 
investigation. 

In  the  course  of  our  investigation  of  spray  combustion,  an  electrical  heating 
system  was  developed  for  the  evaporation  of  the  spray  in  order  to  perform 
experiments  on  the  combustion  of  vapors  of  fuels  such  as  heptane  and  kerosene.  The 
system  could  also  be  used  for  performing  tests  on  the  combustion  of  sprays  in  the 
background  of  their  own  vapors.  Because  of  the  budgetary  constraints,  this  system 
was  only  partially  developed  and  some  preliminary  results  were  obtained  on  heptane- 
vapor  combustion.  A  most  interesting  result  of  these  tests  was  the  occurrence  of  the 
star-shaped  flame  predicted  earlier  by  Sivashinsky  (1983).  Because  of  the 
transparent  quartz  plate  used  in  our  burner  system,  the  geometry  of  the  star  flame 
could  be  clearly  visualized  as  seen  in  the  direct  photograph  shown  in  Fig.  2.  Thus  in 
agreement  with  the  theory,  as  the  rate  of  stretch  was  steadily  increased,  the  cellular 
flame  first  assumed  the  star-shaped  geometry  and  eventually  at  higher  stretch  rates 
the  flame  surface  became  flat.  We  emphasize  ihat  the  geometry  of  the  flame  surface 
relates  to  the  total  flame  surface  area  which  in  turn  controls  the  turbulent  burning 
rates. 

It  is  emphasized  again  here  that  the  results  of  the  basic  research  efforts  outlined 
above  are  for  application  to  modeling  of  turbulent  combustion  process  in  liquid  rocket 
motors.  The  goal  of  these  basic  research  studies  is  to  provide  a  fundamental  insight 
into  the  local  dynamics  of  the  reactive  flow  and  thereby  help  to  improve  our  capability 
for  modeling  the  more  complex  global  process.  For  example,  presently  baffles  are 
being  used  in  the  rocket  motor  cavity  in  order  to  modify  the  standing  acoustic  wave 
patterns  within  the  cavity  and  thereby  prevent  the  coupling  between  chemical  heat 
release  and  the  acoustic  waves.  Studies  that  improve  our  understanding  of  the 
acoustic- wave/flame  interactions  may  be  helpful  in  future  development  of  better 
techniques  for  designing  safe  and  more  economic  rocket  motors.  Also,  it  is 
emphasized  that  the  present  experimental  model  provides  a  stationary  acoustically 
unstable  flame  which  will  facilitate  the  diagnostic  investigation  of  these  unstable 
flames.  Finally,  we  note  that  only  preliminary  data  was  obtained  on  the  acoustic 
frequency  and  the  temperature  fluctuations  of  the  unstable  flames,  Appendix  II. 

Clearly,  more  research  is  needed  to  further  explore  these  and  other  important  aspects 
of  the  unstable  flames  which  is  of  central  significance  to  the  liquid  rocket  combustion 
phenomenon  (Williams,  1985). 


THERMO-DIFFUSIVE  FLAME  INSTABILITIES 

Another  important  type  of  flame-front  instability  in  premixed  systems  is  the 
thermo-diffusive  flame  instability  (Sivashinsky,  1983)  which  is  caused  by  the 
preferential  diffusion  of  heat  versus  reaction-controlling  species.  Thus,  when  the 


Lewis  number  Le  =  a/D  (  a  and  D  refer  to  the  diffusivity  for  heat  and  the  deficient 
component)  of  the  mixture  is  less  than  a  critical  value  Le*  <  1.0,  small  corrugations  of 
the  flame  surface  will  grow,  leading  to  the  formation  of  regularly  spaced  cells  hence 
cellular  or  polyhedral  flames.  Although  the  phenomenon  of  thermo-diffusive  flame 
instability  is  primarily  controlled  by  the  processes  at  the  molecular  scales, 
i.e.molecular  diffusion,  it  is  expected  to  have  a  pronounced  influence  on  the  burning 
rate  in  turbulent  flames  (Williams,  1985).  This  is  in  part  because  of  the  fact  that 
turbulent  flames  are  often  composed  of  an  ensemble  of  stretched  laminar  flamelets 
whose  surface  corrugation  may  be  modified  by  the  thermo-diffusive  instability 
mechanism.  The  influence  of  Le  on  the  turbulent  propagation  velocity  has  been 
discussed  by  Williams  (1985). 

In  a  previous  experimental  study  (Bidinger  and  Sohrab,  1985),  it  was  shown  that 
the  chemical  kinetic  mechanisms  may  play  a  significant  role  in  the  amplification/ 
sustenance  of  thermo-diffusive  flame  instabilities.  This  hypothesis  was  tested  by 
observation  of  the  effects  of  bromo-trifluoro-methane,  CFsBr,  on  polyhedral  butane 
Bunsen  flames.  We  note  here  that  besides  the  early  investigations  by  Markstein 
(1964),  relatively  little  experimental  work  has  been  devoted  to  the  study  of  the 
temperature  profiles  in  cellularly  unstable  flames.  In  particular,  few  experiments  have 
been  performed  on  the  periodic  temperature  variations  in  rotating  polyhedral  flames 
(Sohrab  and  Law,  1985).  Because  of  the  important  fundamental  as  well  as  practical 
value  of  the  problem,  the  influence  of  radical  species  on  cellular  flames  was  re¬ 
examined  and  the  temperature  of  the  trough  and  crest  regions  of  the  polyhedral  flames 
of  butane  were  measured. 

The  results  of  the  experimental  investigation  of  polyhedral  butane  flames  are 
presented  in  the  Appendix  III.  In  accordance  with  the  earlier  observations  (Bidinger 
and  Sohrab  1985),  it  was  established  that  the  radical  species  may  indeed  play  a 
significant  role  in  the  process  of  thermo-diffusive  flame  instabilities.  On  the  basis  of 
the  observations  and  in  terms  of  the  known  radical  scavenging  properties  of  halogens, 
a  phenomenological  mechanism  for  the  observed  effects  of  CT^Br  on  the  polyhedral 
flames  was  presented.  Another  significant  aspect  of  the  study.  Appendix  III,  was  the 
description  of  the  lean  flammability  limit  of  mixtures  on  the  basis  of  the  chemical 
kinetic  mechanisms.  Thus,  it  was  described  how  the  different  temperature 
sensitivities  of  the  important  chain  branching  reaction  such  as 

H  +  02  ->  OH  +  O 

and  a  chain  terminating  reaction  such  as 

H  +  O2  +  M  — >  HO2  +  M 

may  describe  the  occurrence  of  a  critical  minimum  temperature,  and  hence  fuel 
concentration,  below  which  steady  combustion  will  not  be  possible!  Appendix  III). 

The  fundamental  significance  of  the  results  described  above  to  the  understanding 
of  turbulent  combustion  in  liquid  rockets  is  apparent.  It  is  emphasized  here  that  in 
rocket  motors  using  H2  +  O2  as  reactant  ,  the  large  difference  between  the  molecular 
diffusivity  of  the  fuel  versus  oxidizer  makes  the  mixture  particularly  susceptible  to 
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preferential-diffusion  effects.  Almost  every  aspect  of  the  problems  studied  in  the 
Appendix  III  are  in  need  of  further  research.  Chemical  kinetic  description  of  the  limits 
of  flammability  is  of  central  significance  to  the  understanding  of  local  flameiet 
extinction  in  turbulent  premixed  flames.  Also,  it  is  emphasized  that  the  thermo- 
diffusive  instabilities  were  found  to  be  enhanced  in  the  oxygen-rich  environments 
which  are  expected  to  predominate  in  liquid  rocket  motors.  In  turbulent  combustion 
fields,  multiple  flame  surfaces  occur  which  are  in  constant  random  motion  due  to  the 
turbulent  fluctuations  (Williams,  1985a).  Because  of  the  thermo-diffusive  instabilities, 
the  interactions  between  corrugated  flame  surfaces  may  occur  which  will  modify  the 
nature  of  the  reactive  flow  (Sohrab  and  Chao.  1984).  The  central  importance  of 
studies  on  various  mechanisms  of  the  flame  front  instabilities  to  the  understanding  of 
turbulent  combustion  has  been  emphasized  (  Williams,  1985). 


INFLUENCE  OF  VORTICITY  ON  FLAME  STABILITY/STABILIZATION 

Two  hydrodynamic  parameters  which  are  of  central  significance  to  the 
understanding  of  turbulent  combustion  are  the  rate  of  stretch  (Williams,  1985a)  and 
the  vorticity  (Chen,  et  al.,  1987;  Lin  and  Sohrab,  1987;  Sivashinsky  and  Sohrab,  1987). 
In  liquid  rocket  motors,  the  presence  of  two  phases  further  complicates  the  dynamics 
of  the  reactive  field.  Turbulent  eddies  are  expected  to  interact  with  the  flame  surfaces 
as  well  as  the  moving  droplets  in  a  complex  and  transient  fashion  (Faeth,  1983). 
Therefore,the  understanding  of  the  interactions  between  vorticity  and  premixed  and 
diffusion  flames  is  of  central  importance  to  the  global  modeling  of  turbulent  premixed 
and  diffusion  flames.  For  example,  in  recent  theoretical  (Sivashinsky  and  Sohrab, 
1987)  and  experimental  (Chen,  et  al.,  1987;  Lin  and  Sohrab,  1987)  studies,  it  was 
shown  that  flow  rotation  can  reduce  the  minimum  fuel  concentration  at  flame 
extinction.  In  addition,  it  was  found  that  the  rotation  of  the  gas  will  enhance  the 
formation  of  soot  particles  in  opposed-jet  counter-rotating  flows  when  oxvgen- 
enriched  air  is  used.  Such  results  have  direct  relevance  to  the  problem  of  local  flame 
extinction  and  soot  formation  in  turbulent  combustion  environments  of  the  rocket 
motors  using  hydrocarbon  fuel. 

The  influence  of  slow  rotation  on  the  geometry  and  stability  of  flames  propagating 
in  a  rotating  cylindrical  tube  was  theoretically  investigated  in  the  limit  of  weak  thermal 
expansion  (Sivashinsky,  et  al.,  1988).  Thus,  the  effects  of  centrifugal  and  Coriolis 
accelerations  on  the  premixed  flames  were  illuminated.  In  the  study,  a  discrete  set  of 
angular  velocities  were  identified  at  which  resonance  effects  occurred,  leading  to 
unbounded  amplification  of  the  flame  propagation  speed.  Of  course,  in  real  situations 
the  dissipative  viscous  effects  will  prevent  the  occurrence  of  such  resonance  points 
(Sivashinsky  and  Sohrab,  1987).  An  important  prediction  of  the  theory  was 
concerning  the  suppression  of  thermo-diffusive  instabilities,  i.e.  cellular  flame 
structures,  by  the  flow  rotation. which  was  also  substantiated  experimentally 
(Sivashinsky,  et  al.,  1988).  On  the  basis  of  this  result,  we  expect  that  turbulent 
eddies  may  tend  to  suppress  the  flame  corrugations  induced  by  thermo-diffusive 
instabilities.  The  experimental  observation  of  the  shape  of  the  flame  surface  in  'he 
rotating  tube  also  showed  good  qualitative  agreement  with  the  theoretical  predictions. 
The  central  importance  of  these  findings  to  the  understanding  of  turbulent  combustion 


has  been  emphasized  (Sivashinsky,  1983).  Many  other  aspects  of  the  problem  thus 
described  require  further  future  exploration  both  experimentally  as  well  as 
theoretically.  In  particular,  the  theoretical  analysis  may  be  extended  to  include  the 
effects  of  strong  heat  release.  Also,  the  effects  of  preferential  diffusion  and  non-unity 
Lewis  number  require  further  experimental  study. 

In  another  theoretical  investigation  (Sheu,  et  ah,  1989),  the  effects  of  the  flow 
rotation  on  the  geometry  and  stabilization  of  Bunsen  flames  were  investigated.  It  was 
found  that  flow  rotation  first  buckled  the  central  portion  of  the  otherwise  conical  flame 
and  eventually,  at  higher  rotation  velocities,the  stabilization  of  the  flame  on  the  burner 
rim  became  impossible  and  the  flame  flashed  back  into  the  burner  tube  (see  Appendix 
V).  The  predictions  of  the  theory  just  described  were  also  substantiated  by  the 
experimental  observation  of  methane  and  butane  premixed  flames  stabililized  on  a 
rotating  Bunsen  burner.  Both  the  geometry  as  well  as  the  flame  flash-back 
characteristics  observed  experimentally  were  in  good  qualitative  agreement  with  the 
theory.  Many  aspects  of  the  theoretical  and  experimental  investigation  just  described 
are  in  need  of  further  future  exploration.  For  example,  because  of  the  importance  of 
the  flame  stabilization  phenomenon,  the  theoretical  model  may  be  extended  to  include 
the  effects  of  the  heat  loss  to  the  burner  rim.  Also,  additional  experiments  are  needed 
to  improve  our  understanding  of  the  transient  process  of  flame  flash-back.  We 
emphasize  again  that  even  though  the  studies  described  herein  are  necessarily 
confined  to  a  certain  geometrical  model  flow  fields,  their  aim  is  more  general  in  that 
they  may  help  to  improve  our  understanding  of  the  local  flow/flow,  flow/flame  and 
flame/flame  interactions  in  turbulent  reactive  flows. 


NON-PLANAR  FLAME  CONFIGURATIONS 

Flame  stabilization  on  nozzles  is  a  common  occurrence  in  many  combustion 
devices.  In  a  theoretical  investigation  (Sheu  and  Sivashinsky,  1989)  the  non-planar 
solutions  for  nozzle-stabilized  premixed  flames  wc  :e  determined  for  reactive  jets 
impinging  on  a  plane  wall  in  the  stagnation-point  configuration.  The  re-ults  of  the 
theory  (see  Appendix  VI  )were  found  to  be  in  close  qualitative  agreement  with  the 
experimental  observations.  In  particular,  it  was  established  that  in  the  stagnation 
flow  configuration,  besides  the  usual  planar  flame,  other  non-planar  flame  surfaces 
may  also  occur  whose  dynamic  stability  require  further  experimental  and  theoretical 
investigation.  Model  problems  such  as  the  one  being  discussed  here  may  be 
considered  as  the  first  step  towards  the  simulation  of  flame  stabilization  on  injector 
pons  and  flame  holders.  Indeed,  it  is  expected  that  future  stability  analysis  of  the 
cylindrical,  nozzle-stabilized  flames  would  reveal  more  detailed  aspects  of  the  flame 
stand-off  as  well  as  the  flame  blow-off  mechanisms.  More  experimental 
investigations  would  also  be  needed  in  order  to  allow  comparisons  with  the  theoretical 
results  of  the  flame  stability/stabilization  characteristics  just  described. 
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Figure  2  Direct  photograph  or  star  shaped  flame  of  rich  butane-air 
mixture  in  stagnation  -  point  flow  Photographed  from  the 
top  through  quartz  plate(  see  Appendix  T.  Fig  1  ) 
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Abstract — The  steady  combustion  ol'polvdispersed  sprays  of  ethanol  and  kerosene  in  the  stagnation-poini 
rlow  of  lean  methane  air  mixtures  is  experimentally  investigated.  Using  iaser  Doppler  veiocimctry.  the  axial 
and  radial  v  .locity  profiles  of  the  droplets  are  measured  for  water  sprays  in  the  presence  or  absence 
of  lean  methane-air  flame.  These  results  arc  then  compared  with  the  velocity  profiles  obtained  for  small 
MgO  particles  under  identical  flow  conditions.  The  results  show'  the  effects  of  gas  expansion  on  dropiet 
deceleration  t acceleration)  in  the  pre-flame  (post-flame)  regions.  Also,  it  is  found  that  addition  of  water 
sprav  results  in  the  formation  of  a  distributed  region  of  yellow  -orange  emission  dow  nstream  of  the  iean 
methane  flame.  In  combustion  of  cthanoi  spray,  certain  critical  fuci-conccniraiion -velocity  hums  are 
identified  above  which  the  flame  becomes  acousticaiis  unstable.  The  implications  of  tne  stuo\  to  me 
moueiing  of  turbulent  spray  combustion  are  discussed. 

INTRODUCTION 

Studies  on  combustion  of  liquid  fuel  sprays  help  to  improve  the  understanding  of 
burning  processes  in  gas-turbtne.  liquid-rocket  and  diesel  engines,  among  many  other 
important  applications.  The  complications  due  to  the  simultaneous  presence  of  two 
phases  and  uncertainties  in  droplet  size  and  distribution  are  compounded  by  the 
presence  of  turbulence  in  most  practical  applications.  Turbulent  fluctuations  modify 
the  mixing,  evaporation  and  combustion  processes  and  lead  to  transient  interactions 
between  droplet,  flow/flame  dynamics  (Faeth.  1983).  Because  of  the  large  number  of 
parameters,  analytical  and  experimental  studies  often  aim  at  modeling  the  separate 
effects  of  individual  parameters  which  are  known  to  locally  govern  the  reactive 
flow  field.  The  results  of  such  studies  will  help  the  future  development  of  more 
comprehensive  global  me  ds  for  turbulent  spray  combustion. 

The  evolution  of  theoretical  models  of  spray  combustion  has  been  discussed  by 
Williams  (1985).  Recent  developments  in  the  field,  including  the  progress  made 
m  turbulent  sprav  combustion,  have  been  reviewed  by  Faeth  f  1983).  Because  of  the 
comoiexity  of  the  problem,  besides  the  full-scale  studies,  two  distinct  and  compli¬ 
mentary  classes  of  investigations  have  been  found  necessary.  First  are  studies  which 
aim  at  understanding  of  local  phenomena,  and  consider  combustion  of  single  droplet 
or  an  array  of  droplets.  Second,  which  includes  the  work  reported  herein,  are  studies 
concerned  with  the  global  effects  and  investigate  combustion  of  mono-  or  poi>- 
dispersed  sprays  in  laboratory  scale  burners.  For  example,  spherical  propagation  of 
flames  in  mono-dispersed  sprays  of  ethanol  and  u-octane  has  been  studied  in  Wilson- 
cioud-chamber  type  apparatus  (Hayashi  and  Kumagai.  1974:  Hayashi  et  u/..  lq76l.  In 
another  experimental  model,  combustion  of  poly-dispersed  sprays  in  propane- 
kerosene  drop-air  \Mizutani  and  Nakaiima.  1973)  and  kerosene-air  (Polvmeropoulos 
and  Das.  1975)  systems  were  studied  in  inverted-cone-fiame-burner  configuration. 
Also,  the  propagation  of  planar  flames  in  polv-dispersed  sprays  in  cylindrical  vessel 
have  been  studied  (Bailal  and  Lefebvre.  1981).  In  these  studies,  among  other 
important  findings,  the  propagation  speed  of  the  flames  in  sprays  were  determined. 

In  modeling  of  turbulent  spray  combustion,  one  of  the  maior  difficulties  to  be 
addressed  is  the  way  in  which  the  muitipie  flame  surfaces  and  the  individual  droplets 
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behave  under  the  large  velocity  fluctuations  (Faeth,  1983).  Thus,  the  problem  in  part 
involves  interaction  of  a  burning  droplet  with  severe  velocity  non-uniformities.  It  is 
expected  that,  similar  to  homogeneous  turbuient  combustion,  the  dame  stretch 
compression  effects  (Karlovitz  et  a!.,  1953)  will  also  play  an  important  role  in 
turbulent  spray  combustion.  We  note  that  in  passing  through  the  flame  sheet,  droplets 
experience  a  sudden  velocity  jump  induced  by  the  expansion  of  the  gas.  Also,  besides 
the  transient  interactions  between  flamelets  (Sohrab  et  al..  1984;  1986).  mixing  non¬ 
homogeneities  may  lead  to  transitions  between  premixed  and  diffusion  burning 
regimes  t  Lin  and  Sohrab.  1987). 

The  present  experimental  investigation  is  an  extension  of  an  earlier  study  (Sohrab. 
1985)  on  combustion  of  hydrocarbon  sprays  in  laminar  stagnation-point  flow.  The 
study  is  motivated  by  the  fact  that  in  such  flow  configuration,  not  only  the  fluid 
mechanics  of  the  gaseous  flow  is  relatively  well  understood,  but  the  effects  of  stretch 
rate  on  the  flame  can  be  systematically  studied.  Also,  since  steady  combustion  of 
planar  flames  can  be  achieved,  the  investigations  of  the  droplet/flame  interactions  and 
flame  diagnostics  are  facilitated.  In  the  following,  the  stagnation  flow  sprav-burner  is 
described  and  the  axial  and  radial  velocity  profiles  of  the  gas  as  well  as  non-reactive 
dropiets  in  (water  spray-methane-air)  systems  are  presented.  Next,  the  results  on 
combustion  and  acoustic  flame  instabilities  observed  in  (ethanol  spray-methanc-airi 
and  (kerosene  spray-methane-airi  systems  are  discussed.  A  summary  of  our  principal 
findings  is  given  in  the  concluding  remarks. 


SPRAY  BURNER  SYSTEM 

The  stagnation  flow  spray  burner  is  composed  of  a  quartz  contoured  nozzle  with 
33  mm  exit  diameter  and  10  to  1  area  contraction,  a  flat  quartz  plate  and  a  liquid  fuel 
atomizer.  In  Figure  1  the  schematic  of  the  spray  burner  is  shown.  Air  and  gaseous  fuel 
are  metered  by  conventional  rotameters  and  premixed  in  the  low-er  chamber  of  the 
burner  as  shown  in  Figure  I.  The  combustible  mixture  then  enters  the  102  mm 
diameter  extension  tube  containing  small  glass  beads  thus  producing  a  uniform 
velocity  profile  within  the  tube.  A  liquid  fuel  atomizer  is  situated  at  the  center  of  the 
extension  tube  as  shown  in  Figure  1. 

Two  types  of  atomizers  with  fundamentally  different  working  mechanisms  are 
considered  for  generation  of  poiydispersed  sprays.  The  first  type  of  atomizer  employs 
uitrasonic  vibration  of  a  specially  designee  nozzle  using  piezo-eiectric  crystals  for 
liquid  atomization.  The  liquid  fuel  is  fed  into  the  atomizer  by  gravity  from  a  specially 
designed  supply  system  which  maintains  an  optimum  constant  liquid  lex ei  at  ail  times. 
The  atomizer  produces  20-50  urn  droplets  in  a  narrow  jet.  about  4()  degree  cone  angle, 
with  negligible  spillage  and  axial  momentum.  The  second  atomizer  is  a  high-pressure 
atomizer  that  is  used  in  conventional  oil  spray  burners.  The  liquid  is  introduced  to  the 
atomizer  from  a  1  liter  storage  tank  pressurized  to  lOOpsig  by  nitrogen  gas.  The  cone 
angle  of  the  spray  is  about  80  degrees. 

The  droplets  are  convected  by  the  coflowing  methane  air  mixture  as  shown  in 
Figure  1.  In  passing  through  the  contoured  nozzle,  a  uniform-velocity  poiydispersed 
spray  is  formed  at  the  exit  plane  of  the  nozzle.  This  flow  subsequently  impinges 
on  a  flat  quartz  plate  located  at  a  specified  and  fixed  height.  14  mm.  above  the 
nozzle  rim.  As  shown  in  Figure  1.  a  ring  of  cooling  water  with  smail  water  jets 
surrounds  the  nozzle  for  cooling  the  rim  as  well  as  the  exhaust  system.  A  fuel 
drain  pioe  is  provided  for  droplets  which  may  accumulate  within  the  glass  bead 
region.  An  exhaust  system  is  surrounding  the  nozzle.  Figure  I.  to  prevent  tnc 
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FIGURE  ?  Direct  photograph  of  streamlines  in  nonreactivc  sprat  of  water  in  air.  I"  =  :25cm  > 
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FIGURE  -  \u.ii  \eiocil\  proitles  01  water  uropiets  and  MeO  particles  in  reacme  and  non-reacmc  flow > 
!  =  i 25 cm  s.  li  O  ')  Dl^ce  s.  0,  =  =  nmm 


compared  to  particles  (Figure  4t.  On  the  other  hand,  the  radial  velocity  of  droplets 
is  slightly  smaller.  <  2cm  s.  than  that  of  particles  (Figure  5).  Also,  the  small  but  finite 
value  of  the  radial  velocity  aoout  I  5  cm  s.  at  Z  =  0  shows  that  the  flow  at  the  nozzle 
rim  is  not  purely  axial.  This  early  divergence  is  apparcntlx  caused  by  the  overall 
pressure  distribution  in  the  How  held  since  it  remains  nearly  constant  even  in  the 
absence  of  the  exhaust  suction. 
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FIGURE  5  Radial  velocity  protiles  of  water  droplets  and  MgO  particles  in  reactive  and  non-reacme 
flows  I'"  =  125cm  s.  H.O:  O.UIVee  s  Cl,  =  7.7” n .  r  =  bmm. 

The  characteristics  of  water  spray  in  the  presence  of  lean  methane  air  flames 
are  examined  next.  At  nozzle  velocity  of  f  '"  =  124cm. s,  a  lean  methane  air  flame 
with  volumetric  fuel  concentration  Q ,  =  7.5%  is  established  first.  Small  quantity. 
0.02  cc.s.  of  water  spray  is  then  introduced  into  the  flow  field.  The  direct  photograph 
of  the  initial  blue-violet  lean  methane  flame  is  shown  in  Figure  6a.  After  introduction 
of  water  spray,  a  homogeneous  region  of  yellow-orange  radiation  occurs  downstream 
of  the  blue-violet  flame  as  shown  in  Figures  6b-6d.  This  is  quite  interesting,  since  the 
lean  methane  flame  is  not  easily  susceptible  to  soot  formation.  This  phenomenon  is 
in  part  caused  by  the  cooling  effects  of  evaporating  droplets  in  the  hot  post-flame 
regions.  It  is  known  that  the  slow  combustion  of  carbon  monoxide  is  not  complete 
until  about  2  mm  downstream  of  the  flame  zone  (Tsuji  and  Yamaoka.  1982).  Hence, 
the  cooling  effects  due  to  droplet  evaporation  may  suppress  CO  combustion  and  lead 
to  polymerization.  Also,  scavenging  of  radicals,  such  as  OH.  due  to  heterogeneous 
surface  recombination  may  occur.  Hydroxyl  radical  is  known  to  be  effective  in 
combustion  of  soot  precursors. 

Similar  phenomena  were  observed  when  distilled  water  was  used  instead  of  the  tap 
water.  Therefore,  the  yellow-orange  radiation  could  not  be  associated  with  emission 
from  possible  sodium  sait  impurities  in  water.  Small  quantities  of  soot  were  found  to 
deposit  on  thin  metallic  wires  inserted  in  the  yellow  post-flame  regions.  The  yellow- 
orange  radiation  zone  ended  abruptly  about  2  mm  from  the  quartz  plate:  a  position 
corresponding  to  the  edge  of  the  boundary  layer  (Figure  7).  Hence,  no  accumulation 
of  soot  particles  on  the  hot  quartz  plate  was  observed.  The  exact  nature  of  the 
veilow-orange  radiation  region  and  its  possible  relation  to  soot  formation  requires 
further  future  exploration.  Figures  6e-6d  are  direct  photographs  of  droplet  stream¬ 
lines  that  are  illuminated  by  a  narrow  vertical  laser  sheet-light,  respectively  upstream 
and  downstream  of  the  flame  zone.  Comparisons  of  Figures  6a  with  6b-d  show  that 
the  flame  thickness  increases  in  the  presence  of  water  spray. 
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FIGURE  6  Dircci  photograph  of  CH,  air  flames  in  presence  of  scaler  droplets  la)  Pure  CH.  air  flame 
lb)  CH,  air  flame  uuh  water  drops  ic)  CH.  air  flames  with  scaler  drops  under  pre-flame  illumination  (d) 
CH,  air  rlames  ssim  ssater  drops  and  post-flame  tilunnnaiion  I  =  124cm  s.  £2.  =  '  NV  H  .O  0.02  ec  - 

The  axial  temperature  profiles  for  the  lean  methane  air  llatnes.  Q,  =  T."“»  and 
f  =  1 25 cm  .s.  with  and  without  water  spray  are  shown  m  Figure  The  data  are 
obtained  by  axial  traverse  oi  a  siiica-coated  Pt-Rt-i7"„  Rh  thermocouple  witn 
0.1  mm  wire  diameter  at  a  constant  radial  position  r  —  nmm.  Cieariy.  the  measured 
temperatures  represent  a  mean  value  since  droplet  collisions  with  the  thermocouple 
lunction  occur.  The  cooling  elTects  of  water  spray  reduce  the  tlame  propagation  speed, 
which  is  manifested  r>\  the  relocation  ol  tlame  closer  to  the  stagnation  piane. 
Z  =  14  mm.  The  axial  thickness  of  the  orange-yellow  layer  is  identified  in  Figure  “. 
The  upper  surface  of  the  yellow  emission  layer  abruptly  ends  at  a  finite  distance  from 
the  hot  quartz  plate  (Figure  6b i.  it  is  noted  that  relatively  small  changes  of  the  mean 
temperatures.  <  75°C.  have  resulted  in  substantial  change  in  radiation  characteristics 
in  post-flame  regions.  Of  course,  because  of  the  vaporizing  droplets,  the  local  tem¬ 
perature  fluctuations  are  much  larger  than  75°C. 

The  variations  of  the  axial  and  radial  velocity  as  functions  of  Z  lor  droplets  and 
MgO  particles  in  the  presence  of  lean  methane  air  tlame  are  also  shown  in  Figures  4 
and  5.  As  seen  in  Figure  4.  the  differences  between  axial  velocity  of  particles  and 
droplets  are  quite  substantial.  The  data  show  that  droplet  deceleration  (acceleration i 
induced  by  the  gas  expansion  m  the  pre  (post)  flame  zones  is  smaller  than  that 
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FIGURE  7  Axial  temperature  profiles  in  lean  methane.air  flames  (a)  without,  (b)  with  water  spra> 
Q,  =  7.7%.  r"  =  I25cm/s.  H.O  =  0.019ccs.  r  =  6mm. 

corresponding  to  particles.  While  the  particle  velocity  V  jumps  from  30  to  120cm  s 
across  the  flame  front,  that  of  droplets  is  from  40  to  80  cm/s.  This  is  to  be  expected 
in  view  of  the  larger  inertia  of  the  droplets  and  their  higher  resistance  to  velocity 
changes  as  compared  to  the  particles.  That  is.  since  the  droplets  are  heavier  than 
the  MgO  particles,  their  acceleration  by  the  expanding  gas  will  be  more  difficult. 
Obviously,  such  velocity  jumps  result  in  enhanced  droplet  slip  velocity  which  will 
substantially  alter  droplet  evaporation/combustion  in  chemically  reactive  sprays. 

In  the  stagnation-point  flow  configuration,  the  axial  velocity  at  the  upstream  edge 
of  the  preheat  zone  may  be  defined  as  the  laminar  propagation  speed  of  the  stretched 
flame  which  differs  from  the  propagation  speed  in  the  absence  of  stretching.  Accord¬ 
ing  to  Figure  4.  this  axiai  velocity  is  about  30 cm/s  at  the  position  Z.  ^  S  mm.  based 
on  the  MgO  velocity  profile,  in  agreement  with  prior  observations  for  lean  methane 
flames  (Yamaoka  and  Tsuii.  1984:  Wu  and  Law.  1984)  under  similar  stretch  and 
concentration  conditions.  Because  of  the  inertial  effects,  the  velocity  profile  based  on 
droplet  motion  differs  from  that  of  the  gas  and  cannot  be  used  to  deduce  the  flame 
speeds  in  the  presence  of  water  spray.  However,  for  the  flame  cooled  by  water  spray, 
an  approximate  and  lower  propagation  speed,  about  20crms.  is  estimated  when  the 
velocity  profile  based  on  MgO  particles  is  linearly  extrapolated  to  the  corresponding 
flame  position  (Figure  4). 

The  differences  between  axial  distribution  of  radial  velocity  for  particles  versus 
droplets,  shown  in  Figure  5.  also  reflect  the  effects  of  higher  droplet  inertia.  1 1  is  noted 
that  the  boundary  layer  thickness  d  for  reactive  flow  is  about  twice.  b  =  2  mm.  that 
of  the  cold  flow  (Figure  5).  This  tendency  is  anticipated  since  for  the  present  flow  field 
b  is  proportional  to  v'v  (Batchelor.  1970)  and  the  kinematic  viscosity  v  is  an  almost 
linearly  increasing  function  of  temperature.  In  addition,  because  of  the  composition 
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differences,  the  value  of  v  is  different  for  the  reactive  and  non-reacuve  flows.  The  data 
in  Figure  5  also  show  that  the  measured  maximum  radial  velocities  of  droplets  and 
particles  are  nearly  the  same.  66  and  68  cm/s,  while  those  within  the  boundary  layer 
are  identical. 


REACTIVE  SPRAY  IN  REACTIVE  FLOW 

In  this  section,  the  combustion  of  fuel  sprays  in  lean  methane-air  mixtures  are 
discussed.  Basically,  two  distinguishable  burning  modes  may  be  identified  depending 
on  the  volatility  of  the  liquid  fuels.  For  less  volatile  fuels,  such  as  kerosene,  droplets 
remain  intact  until  they  closely  approach  the  methane  flame.  Hence,  evaporation 
occurs  primarily  in  the  preheat  zone  and  downstream  of  the  reaction  zone.  For  more 
volatile  fuels,  such  as  heptane  and  ethanol,  on  the  other  hand,  droplet  evaporation 
begins  immediately  after  atomization.  Hence,  in  this  case  smaller  droplets  burn  in 
the  gaseous  mixture  of  their  own  vapor  and  methane-air.  For  heptane  sprays,  the 
concentration  of  methane  in  air  could  be  reduced  to  zero  while  maintaining  steady 
combustion  of  the  spray. 

Direct  photographs  of  flames  in  sprays  of  ethanol  and  kerosene  have  been  obtained 
and  typical  samples  are  shown  in  Figures  8  and  9.  respectively.  In  Figure  8.  the  side 
view  as  well  as  the  top  view,  through  the  quartz  plate,  are  shown.  A  horizontal  sheet 
of  laser  light  was  axially  positioned  5  mm  upstream  of  the  flame  surface  in  Figure  Sb 
to  illuminate  the  droplets.  It  is  evident  that  with  higher  intensity  and  narrower 
thickness  of  the  laser  sheet-light,  such  photographs  will  allow  for  approximate  evalu¬ 
ation  of  droplet  size  and  distribution  in  the  pre  post-flame  regions.  Preliminary 
observations  of  magnified  photographs,  such  as  Figure  8b.  showed  that  most  of  the 
droplets  are  not  spherical  in  shape.  For  kerosene  sprays.  Figure  9.  envelope  flames  are 
observed  to  surround  the  individual  droplets  downstream  of  the  flame  front.  Also, 
severe  sooting  and  eventual  accumulation  of  soot  on  the  quartz  plate  occurs  as  the 
kerosene  flow  rate  is  increased.  The  photographs  for  kerosene  sprays  were  taken  when 
the  high-pressure  atomizer  was  being  used  (Sohrab.  1985).  It  was  found  that  with 
dense  sprays,  kerosene  drops  actually  impinge  on  the  hot  plate  and  support  a  partially 
premixed  flame  (Sohrab  et  al..  1984)  adjacent  to  the  hot  piate.  In  this  burning  regime, 
a  thick  layer  of  soot  develops  on  the  quartz  plate  which  is  difficult  to  remove. 

The  variation  of  axial  flame  position  Z.  and  the  mean  flame  temperature  T.  with 
methane  concentration  in  airfL  forethanoi  and  kerosene  sprays  are  shown  in  Figures 
10  and  1 1.  In  these  tests,  the  total  liquid  flow  rate  into  the  burner  is  held  tixed  by 
maintaining  a  constant  atomizer  frequency  and  liquid  level  in  the  supply  tans. 
However,  the  quantities  of  liquid  fuels  which  actually  reach  the  flames  will  depend  on 
the  gas  velocity  I'".  The  effective  values  of  liquid  flow  rates  given  in  Figures  10  and 
11.  are  obtained  by  subtraction  of  the  measured  rates  of  liquid  drainage  from  the 
burner  from  those  introduced  into  the  atomizer.  Such  experiments  are  difficult  since 
the  system  must  reach  a  steady  stale  operation  before  the  measurements  are  made. 
Moreover,  with  the  ultra-sonic  atomizer,  our  ability  to  prescribe  a  certain  desired  rate 
of  liquid  flow  is  limited.  Such  limitations  could  be  overcome  by  the  use  of  more 
advanced  aerosol  generators. 

For  both  ethanol  and  kerosene  sprays,  the  flames  extinguished  below  a  minimum 
methane  concentration  Q,  <  2%.  Hence,  under  the  present  rate  of  stretch  and  droplet 
size.  20-50 pm.  steady  flame  propagation  in  such  sprays  is  not  possible.  The  luminous 
flame  zone  in  the  presence  of  sprays  appeared  thicker  than  that  of  homogeneous  lean 
methane, -air  flame,  about  I  mm.  The  variation  of  Z,  with  flL.  Figures  10  and  I  !.  is  in 
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FIGURE  S  Direct  photographs  of  lean  CH4  air  flames  with  ethanol  spray  ta)  side  view  thi  top  view  »nh 
.ascr  Cncct  iiuhtinc. 


accordance  with  the  known  dependence  of  flame  propagation  speed  on  the  :m\ture 
stoichiometry.  For  ethanol  spray.  Figure  !0.  at  -  3.77  the  flame  flashes  back  into 
the  nozzle.  We  note,  however,  that  at  l'“  =  121  cm  s.  flame  flash-back  does  not  occur 
for  methane  air  flames  even  under  stoichiometric.  Q,  =  9.5.  conditions.  Moreover, 
the  flame  propagation  speeds  in  the  range  45-25  cm  s  have  been  reported  for  pure 
ethanol  sprays  (Hayashi  and  Kumagai.  1974:  Havashi  ei  u/..  1976).  Therefore, 
the  observed  phenomenon  shows  that  higher  laminar  flame  propagation  speed  has 
occurred  in  the  heterogeneous  system  as  compared  to  those  achievable  in  homogen¬ 
eous  mixture  of  either  fuel.  Such  tendencies  have  been  observed  earlier  (Mizutam 
and  Nakaiima.  1973a:  Polvmeropoulos  and  Das.  1975)  when,  under  fixed  overall 
stoichiometry,  the  flame  propagation  speed  was  found  to  first  increase  and  then 
decrease  with  droplet  diameter.  This  tendency  is  in  agreement  with  the  early  pre¬ 
dictions  by  Williams  (1960). 

Unfortunately,  the  variation  of  flame  propagation  speed  with  Q.  in  sprays  could 
not  be  measured  because  of  the  difficulty  in  obtaining  the  axial  velocity  of  the  gas  in 
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FIGURE  9  Direct  photograph  of  lean  CHj  air  flame  with  kerosene  spray 

reactive  sprays.  Here,  the  laser  veloctmetry  based  on  droplet  motion  cannot  be  us-*d 
since  droplets  do  not  follow  the  gas  (Figure  4).  On  the  other  hand,  addition  of  MgO 
particles  to  the  sprays  is  not  feasible,  since  particles  tend  to  agglomerate  within  the 
liuuid  fuel.  When  values  of  Z.  in  Figures  10  and  1 1  are  linearly  extrapolated  to  id,  —  o 
limit,  it  is  found  that  Z.  is  about  14  mm.  ..c.  position  of  (he  stagnation  piano.  This 


FIGURE  iO  Flame  temperature  and  position  js  a  I'uncnon  of  fl,  for  ethanoi  spray  :n  lean  methane  air 
mixtures.  Ethanol  flow  rate  =  0.062  cc  s. 


-  •  4 


FIGURE  1 1  Flame  temperature  and  position  as  function  ol  Q,  for  kerosene  sprat  in  lean  methane  air 
mixtures  Kerosene  flow  rate  =  OOl35cc  v 


tendency  suggests  that  the  flame  propagation  speed  in  pure  sprays,  in  absence  of 
methane,  is  zero  which  corresponds  to  extinction  condition,  hence,  in  absence  of 
sufficient  gas-phase  fuel,  the  gasification  rate  in  spray  is  unable  to  support  self- 
sustaining  flames  in  ethanol  and  kerosene  sprays  considered  herein. 
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An  interesting  observation  made  :n  connection  with  the  tests  on  reactive  -prays  wa- 
the  occurrence  of  acoustically  unstable  flames  in  certain  range  of  nozzle  velocity  and 
fuel  concentration.  For  example  with  1'"  =  174cm  s  and  ethanol  flow  rate  of 
0.08 5cc,s.  it  was  found  that  near  ClF  4%.  the  dame  started  to  vibrate  and  emit 
audible  acoustic  sound.  Such  acoustically  unstable  flames  persisted  beyond  QF  -  4°  o 
until  the  eventual  extinction  of  flame  at  the  rich  fiammibility  limit.  The  acoi  siic  flames 
showed  small.  <2  mm.  axial  oscillations  at  high  frequency,  about  120  Hz.  measured 
from  amplified  microphone  signal.  The  extent  to  which  hydrodynamic,  thermo- 
diffusive  and  chemical  flame  instability  mechanisms  are  effective  in  promotion  and 
sustenance  of  such  acoustic  flames  are  as  yet  not  understood.  The  hydrodynamic  and 
thermo-diffusive  aspects  of  premixed  flame  instabilities  have  been  discussed  bv 
Sivashinsky  (1976).  It  appears  that  resonance  effects  with  the  acoustic  waves  within 
the  burner  cavity  play  a  dominant  role  in  the  present  situation  (Zinn.  I9S6).  Similar 
acoustically  unstable  flames  were  also  observed  earlier  in  kerosene,  heptane,  octane 
and  methanol  sprays  (Sohrab.  1985).  Since  the  onset  of  acoustic  instability  was  found 
to  be  repeatable,  the  present  experimental  model  may  be  used  for  systematic  study  of 
this  interesting  phenomenon. 


DISCUSSIONS 

In  the  sequel,  the  relevance  of  some  of  the  results  obtained  for  non-reacme  and 
reactive  sprays  to  the  global  aspects  of  turbulent  spray  combustion  will  be  described. 
As  in  many  applications,  the  spray  is  assumed  to  be  injected  into  a  hot  and  turbulent 
oxidizing  atmosphere.  Thus,  from  the  central  core  of  the  spray  to  the  far  field  regions 
of  the  atmosphere,  an  entire  spectrum  of  compositions,  from  pure  fuel  to  pure 
oxidizer  are  expected  to  occur.  Since  chemical  reactions  are  considered  to  occur  only 
in  the  gas  phase,  in  the  following  discussions  the  concentrations  refer  to  those  in  the 
gaseous  phase.  We  may  then  expect  that  in  addition  to  regions  of  pure  fuel  <F)  and 
oxidizer  (O).  other  intermediate  premixed  regions  composed  primarily  of  fuel 
(F  —  rO)  or  oxidizer  (O  +  <:F)  will  also  occur.  Finally,  there  will  be  regions  where 
near  stoichiometric  mixtures  (F  —  O)  are  encountered.  Such  regions  may  contain  fuel 
droplets  and  will  be  under  constant  random  motion  due  to  turbulent  fluctuations. 
Hence,  the  composition  of  various  regions  may  change  due  to  droplet  evaporation 
and  by  fast  turbulent  mixing  or  the  relatively  slow  molecular  diffusion.  However, 
ander  stiieiiy  homogeneous  turbulent  situations,  the  stochastic  averages  in  the  multi¬ 
phase  held  may  remain  invariant  in  time. 

In  the  presence  of  chemical  reactions,  some  of  i he  interfaces  between  the  region- 
identified  above  may  support  diffusion  flames  (DF)  or  partially  premixed  1PPF1 
flames  (Sohrab  et  a/..  1984.  1986).  Also,  in  the  regions  ( F  4-  Ol.  premixed  flames  ( PF) 
may  occur  which  will  be  dynamic  due  to  their  propagation.  The  diffusion  flames, 
w  hich  do  not  propagate,  may  be  convected  by  the  motion  of  burning  droplets  or  by 
the  bulk  motion  of  the  enure  interface  between  (Fi  and  (O)  regions.  The  existence  of 
these  multiple  flame  surfaces  is  known  to  depend  on  both  local  fuel  and  oxidizer 
concentrations  as  well  as  the  local  rate  of  stretch.  The  critical  limits  of  fuel  and 
oxidizer  concentration  for  flammability  of  (DF)  and  (PF)  have  been  extensively 
studied.  Moreover,  the  influence  of  the  rate  of  stretch  on  the  flame  extinction  has  been 
discussed  (Williams.  1985).  Thus,  the  combustion  field  may  be  composed  of  many 
(PF).  (DF)  and  (PPF)  surfaces  that  are  in  random  motion  as  droplets  cross  such  flame 
surfaces.  Clearly,  to  understand  the  complex  dynamics  of  the  flow  flame  droplet 
interactions,  simple  models  of  the  local  flow  conditions  should  be  studie  1  first. 
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In  view  ot'ihe  above  considerations,  it  is  evident  that  the  study  of  spray  combustion 
in  stagnation  How  will  help  to  simulate  the  local  How  conditions  within  the  com¬ 
bustion  field.  For  example,  the  problem  of  non-rcactive  spray  discussed  earlier  is  also 
relevant  to  situations  where  fuei  droplets  pass  through  rich  premixed  flame  surfaces. 
Under  such  conditions,  only  droplet  evaporation  and  the  consequent  cooling  effects 
will  occur,  since  downstream  of  the  rich  (PF)  no  more  oxygen  will  be  available  for 
further  combustion.  We  note  however  that  the  behaviour  of  fuel  droplets  after 
rich  (PF)  will  be  somewhat  different  than  water  droplets  since  the  former  may 
undergo  pyrolysis.  Nonetheless,  similar  to  the  observed  effects  of  w  ater  sprays  on  lean 
methane-air  flames,  fuel  sprays  in  fuel  rich  mixtures  are  also  expected  to  influence  the 
process  of  soot  formation  in  the  post  flame  regions. 

The  knowledge  of  velocity  and  trajectory  of  droplets  as  they  approach  and  pass 
through  flames  is  quite  important  for  the  understanding  of  turbulent  spray  com¬ 
bustion.  For  example,  the  observed  enhancement  of  droplet  slip  velocities  across 
flame  fronts  are  expected  to  result  in  local  strained,  stagnation  flow  fields  in  the 
proximity  of  the  individual  droplets.  Such  local  stretch  rates  will  then  alter  the 
structure  and  extinction  behavior  of  the  diffusion  flames  (Williams.  !9"5)  which  mu> 
surround  the  droplets.  Another  important  problem  to  be  addressed  is  the  evaluation 
of  the  critical  minimum  concentration  of  gaseous  fuel  which  would  insure  steady 
combustion  within  the  spray.  For  example,  the  combustion  of  (ethanol  spray-methane- 
air)  system  discussed  earlier  was  symbiotic  in  that  individually,  both  t methane-air) 
and  (ethanol-spray-air)  mixtures  were  incapable  of  supporting  self-sustained  flames 
under  the  specified  flow  conditions.  The  relevance  of  studies  on  spray  combustion  in 
counterflows  to  turbulent  spray  combustion  modeling  were  discussed  in  more  detail 
elsewhere  (Sohrab.  I9S5/. 

CONCLUDING  REMARKS 

In  this  study,  the  behavior  of  reactive  and  non-reactive  sprays  within  methane  air 
mixtures  in  the  stagnation-point  flow  configuration  were  investigated.  In  particular, 
the  motion  of  water  sprays  in  lean  methane  air  flame  was  studied.  It  was  found  that 
water  spray  resulted  in  the  formation  of  a  yellow-orange  radiation  layer  downstream 
of  the  flame  zone.  The  axial  and  radial  velocity  profiles  of  water  droplet  and  small 
MaO  particles  were  measured  in  both  reactive  and  non-reactixe  flows.  Thus,  the  effect 
of  gas  expansion  on  the  droplet  motion  in  the  pre-  and  post-flame  regions  were 
determined  and  compared  to  those  for  MgO  particles. 

For  combustion  of  ethanol  and  kerosene  sprays  in  methane  air  mixtures,  the 
variation  of  flame  temperature  and  position  with  methane  concentration  was  deter¬ 
mined.  It  was  found  that  beiow  a  minimum  methane  concentration,  steady  flame 
propagation  in  the  sprays  was  not  possible.  Also,  above  certain  critical  methane 
concentrations,  the  flames  became  acoustically  unstable.  This  interesting  and  complex 
phenomena  relates  to  acousuc-chemically  coupled  flame  instabilities  and  requires 
further  systematic  investigation.  Since  steady  and  stationary  flames  were  achieved,  the 
present  experimental  model  provides  a  good  vehicle  for  future  investigations  of  spray 
combustion. 
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INTRODUCTION 

In  almost  all  practical  applications,  such  as  liquid  rocket  mocors, 
diesel  engines  and  gas-curbines ,  spray  combustion  occurs  in  curbulenc  flows. 
3ecause  of  che  large  number  of  parameters,  and  interactions  between 
flow/flame/droplet,  the  discripcion  of  the  combustion  field  is  exceedingly 
complex  [1,2].  Therefore,  most  modelling  techniques  first  actempc  to  identify 
the  role  of  che  individual  parameters  which  influence  the  combustion  field  and 
investigate  their  separate  effects  in  simple  analytical  models  or  controlled 
laboratory  experiments.  Such  information  may  then  be  collectively  applied  to 
che  development  of  comprehensive  models  with  proper  measures  taken  to  include- 
che  Interactions  between  the  individual  processes. 

Two  important  hydrodynamic  parameters,  among  others,  which  are  known  to 
influence  the  turbulent  combustion  are  the  rate' of  stretch  and  vorticity.  The 
importance  of  flame  stretch  in  turbulent  combustion  within  homogeneous  gaseous 
flows  is  well  established  [2].  Relatively  less  is  known  about  che  role  of 
vorticity  in  flame-front  dyanmics  [3-5].  In  turbulent  sprays,  on  the  other 
hand,  che  importance  of  the  interaction  between  droplets  and  turbulent  eddies 
has  been  recognized  [1],  In  general,  turbulent  fluctuations  will  induce  flame 
stretching  and  effect  che  motion  of  droplets.  For  example,  the  nature  of 
droplet  evaporation  and  drag  force  on  drolets  are  known  to  depend  on  relacive 
speed  between  che  droplet  and  gas  {  1  ] .  Hence,  studies  on  che  effects  of 
stretch  on  spray  combustion  in  laminar  flows  may  help  the  understanding  of  che 
mechanisms  through  which  both  che  flame  and  droplets  improvise  for  large 
velocity  gradients  which  predominate  che  curbulenc  sprays. 

In  recent  experimental  investigations  [6,7],  the  combustion  of  hydro¬ 
carbon  sprays  in  stagnation-point  flow  was  studied.  Such  flow  configuration 
is  best  suited  for  systematic  study  of  the  effects  of  controlled  velocity 
nonunif ormity  on  flames  [2]  .  In  application  to  spray  combustion,  since  steady 
and  planar  flames  can  be  stabilized  within  the  flow,  che  flame  diagnostics  and 
evaluation  of  droplet  trajectory  through  the  flame  are  facilitated.  In  che 
presenc  experimental  study,  certain  combustion  characteristics  of  ethyl- 
alcohol  3prays  in  lean  mechane/air  mixtures  are  investigated. 

EXPERIMENTS  AND  DISCUSSIONS 

The  3pray  combustion  is  studied  in  a  burner  which  is  schematically  shown 
in  Fig.  1  and  was  described  previously  [6,7].  The  polydispersed  spray  of 
ethyl-alcohol,  20-50  u  droplets,  Is  generated  by  an  ultra-sonic  acomizer  which 
is  surrounded  by  uniform-velocity  flow  of  lean  methane-air  mixture.  After 
passage  through  a  contoured  nozzle,  the  spray  impinges  on  a  planar  quartz 
piace.  The  trajectories,  flow  streamlines,  of  the  iropiets  are  scnematically 
clllustraced  in  Fig.  I.  .Also,  direct  phocograon  of  the  streamlines  nave  oeen 


made  using  laser  sheet-lighting  for  ilumination.  Planar  f lanes  can  be 
obtained  in  Che  echyl-alcohol  spray  with  small  volumetric  concentration  2  of 
methane  in  air,  *  2. 

A  constant  Cocal  flow  race  or  acomizea  liquid,  Q  *  0.0833  cm^/sec.,  is 
maintained  by  controlling  the  level  of  fuel-supply  and  the  ultra-sonic  fre¬ 
quency  of  the  atomizer.  Air  and  methane  flow  rates  are  controlled  by  conven¬ 
tional  rocamecers.  Using  a  pc-pt-lGZ  Rh  thermocouple,  .10  am  wire,  the  tem¬ 
perature  T«  and  position  Ze ,  measured  from  nozzle  exit  plane,  of  the  planar 
flame  are  then  measured  as*2  is  increased  at  fixed  air  flow  rate.  Hence,  the 
mean  gas  velocity  ac  the  exit  plane  of  the  nozzle  V  is  nearly  constant  at 
157  ca/s.  In  this  flow  configuration,  the  flame  poisiton  can  be  related  to 
the  laminar  flame  propagation  velocity  Se  for  stretched  flames  [7].  The 
results  of  such  measurements  for  pure  methane/air  flames  and  those  in  presence 
of  echyl-alcohol  spray  are  shown  in  rigs.  3  and  4.  As  is  to  be  expected,  Zc 
decreases  (Sc  increases)  with  2  until  the  mixture  has  reached  an  affective 
stoichiomecric  composition,  2 ,  *  9.5  and  5.4  in  Figs.  3  and  ■*,  and  tnereafter 
increases  (decreases)  for  fuel-rich  mixtures. 

An  interesting  observation  made  in  the  tesc  was  concerning  the  occurence 
of  certain  acoustic  instabilities  in  both  pure  lean  mechane/air  flames  and 

those  with  spray  of  ethyl-alcohol.  It  is  found  chat  the  planar  flames  start 

vibrating,  emitting  sound  waves,  over  certain  composition  range,  2^  <  2  <  , 

identified  in  Figs.  3  and  4.  Some  of  the  acoustic  frequencies  as  well  as  the 

vibration  frequencies  of  the  flames  have  been  measured.  More  detailed  study 
of  the  complex  coupling  between  the  hydrodynamic  flow  instabilities  and  the 
incrisic  flame  instabilities  [9]  are  presently  being  further  investigated. 
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A3 S TRACT 


The  temperature  fields  of  stationary  and  rotating  polyhedral  3unsea 
flames  of  butane-air  are  experimentally  investigated.  As  the  additives  C7->3r 
and  N«  are  gradually  introduced  into  the  combustible  mixture,  the  variation  of 
the  temperature  in  the  trough  and  crest  regions  of  the  polyhedral  flames  are 
examined.  It  is  found  chat  addition  of  0.36  (5.0)  molar  percentage  of  CF~3r 
(N2)  completely  removes  the  cellular  structure  of  the  polyhedral  flame.  The 
results  suggest  thac  in  conjunction  with  the  preferential-diffusion  mechanism, 
chemical-kinetic  effects  also  play  some  role  in  the  formation,  amplification 
or  sustenance  of  cellular  flame  structures.  Some  phenomenological  arguments 
which  are  based  on  kinetics  of  branched-chain  reactions  are  presented  in 
support  of  the  emoirical  observations.  .Also,  in  butane/air  mixtures  of 
different  velocity,  the  critical  flame  temperatures  corresponding  to  the  onset 
of  ceil  formation  are  found  to  be  nearly  identical.  The  periodic 
circumferential-fluctuations  of  the  temperature  at  different  radial  positions 
within  rotating  polyhedral  flame  of  butane  are  measured  and  the  angular  speed 
of  rotation  is  determined. 
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INTRODUCTION 


Investigation  of  flame-front  instabilities  has  been  recognized  as  a 
significant  branch  of  the  combustion  science  (Williams,  1985).  In  the  purely 
fluid  mechanical  systems,  the  flow  stability  is  related  to  the  Reynolds  number 
which  describes  the  eventual  onset  of  turbulence.  If,  In  addition  to  the 
velocity  gradients,  both  temperature  and  concentration  gradients  are  also 
present,  the  flow  stability  is  further  complicated  due  to  the  simultaneous 
transport  of  momentum,  heac  and  mass.  Thus,  inclusion  of  two  additional 
parameters,  namely  Prandti,  ?r  (ratio  of  momentum  to  heac  diffusivicy)  and 
Lewis,  Le  (ratio  of  heac  to  mass  diffusivicy)  numbers  will  be  needed. 

Finally,  if  the  non-homogeneous  flow  field  is  also  chemically-reaccive,  the 
chemical-kinetic  effects  also  induce  their  direct  or  indirect  influences  on 
the  flow/ flame  stability.  Thus,  in  general  the  stability  of  flames  may  be 
said  to  be  governed  by  hydrodynamic,  thermal,  dixfusional  as  well  as  chemical 
effects.  One  of  the  objectives  of  the  present  investigation  is  to  ascertain 
the  extent  to  which  the  last  category,  namely  chemical-kinecic  effects, 
influences  the  stability  of  premixed  Bunsen  flames,  such  as  polyhedral  flames. 

The  first  report  about  polyhedral  flames  (?HF)  was  given  by  Smirheils  and 
Ingle  (1892).  They  reported  the  existence  of  stationary  or  rotating 
polyhedral  flames  in  fuel-rich  mixtures  of  higher  hydrocarbons  such  as 
benzene,  pentane  and  heptane  burning  in  air.  Lacer,  Smith  and  Pickering 
(1929)  conducted  systematic  experiments  on  ?3F  in  rich  propane/ air  and 
observed  3  to  7-sided  polyhedrals  and  determined  the  characteristic  cell  size 
of  the  flames.  Subsequently,  harkstein  (1964)  performed  comprehensive  studies 
on  the  stability  and  structure  of  cellular  and  polyhedral  flames.  Also,  the 
concentrations  of  soma  of  the  stable  species  in  the  trough  and  crest  zones  of 
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?HF  have  been  measured  'Jose,  ec  al.,  1932;  Marksceir. ,  152  3).  .'.ora  recently. 
Che  influence  of  Che  bumer-rin  aerodynamics  on  PHT  of  bucane/air  on  Bur.sen 
burner  was  scudied  (Sohrab  and  law,  1985)  where  some  preliminary  daca  on  me 
angular  velocity  of  rocating  ?HF  was  also  presented. 

Mancon,  ec  al.  (1952)  attributed  the  occurrence  of  ceiluiar  structures  to 
the  preferential  diffusion  of  'neat  versus  the  deficient  reactant.  Thus,  for 
La  (based  on  the  diffusion  coefficient  of  the  deficient  component)  less  than 
unity,  the  fiame-fronc  corrugations  would  grow  because  of  the  favorable 
stratification  of  the  combustible  mixture  (Lewis  and  von  Elbe,  1987). 

The  theoretical  investigations  of  ceiluiar  flames  were  initiated  by  the 
introduction  of  the  dif f usionai-chern- 1  stability  node!  (3arenblacc,  ec  al., 
1963;  Sivasninsky ,  1976)  and  have  shown  general  qualitative  agreement  with  the 
experimental  observations  (haricstain,  I96«i).  Recently,  more  comprehensive 
models  which  include  Che  influences  of  both  hydrodynamic  as  well  as  thermo- 
diffusive  effects  have  been  reported  (Ciavin  and  Williams,  1982;  Macaion  and 
Matkowsky,  1982).  The  linear  stability  of  non-adiabatic  flames  (Jouiin  and 
Ciavin,  1979;  Sohrab  and  Chao,  1981)  3nd  flames  with  cwo-reactanc  chemistry 
modei  (Jouiin  and  hitani,  1982)  have  also  been  analyzed. 

Although  nose  of  the  observations  of  cellular  flames  thus  far  do  agree 
•with  the  general  hypothesis  of  the  preferential-diffusion  mechanism,  there 
remain  a  few  discrepancies  which  warrant  further  considerations.  First, 
ceiluiar  structures  were  observed  in  rich  methane/ air  flames  (3ocha  and 
Spalding,  1959)  which  is  noc  in  accordance  with  the  preferential-diffusion 
axioms.  On  the  ocher  hand,  it  was  reported  elsewhere  (3ehrens,  1953)  that  to 
polyhedral  flames  could  be  obtained  in  aechane/air  mixtures.  .Also,  cellular 
structures  have  been  observed  in  diffusion  flames  (Garsiae  and  Jackson,  1953) 
which  do  not  possess  any  propagation  velocity  and  as  such  do  not  fail  into  the 
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framework  of  che  model  developed  for  or  as  axed  ays  cans .  Fir.allv,  cr.e  pnvsical 
extant  and  magnitude  of  Che  observed  pnenomena  appear  to  be  coo  large  co  be 
accounted  on  che  basis  or  the  snail  deviation  of  Le  from  unity. 

In  view  of  che  above  consideration,  che  role  of  chemical-kinetic  effects 
through  che  active  radical  species  in  the  formation,  ampiif icacion  and 
suscencanca  of  flame-front  corrugations  was  recently  studied  by  Bidinger  and 
Bohra'o  (1986).  The  results  of  this  study  showed  that  radical  species,  such  as 
atomic  hydrogen,  indeed  piay  some  role  in  che  thermo-diffusive  flame 
instabilities.  Therefore,  it  appears  chac  a  complete  description  of  the 
flame-front  stability  may  require  the  inclusion  of  chemical  effects  in 
addition  to  the  hydrodynamic,  thermal  and  diffusional  effects.  The  present 
investigation  is  an  extension  of  the  previous  study  (Bidinger  and  Bohra'o, 

1986)  and  aims  to  further  explore  the  chemical-kinetic  and  thermo-diffusive 
nature  of  cellular  flame  structures.  Such  effects  are  explored  by  examination 
of  temperature  in  che  trough/ crest  cones  when  either  chemical  inhibitor  Cf-.Br 
or  nitrogen  are  gradually  added  to  PHF  of  rich  butane/air  3unsen  flames. 

Also,  the  rotation  speed  and  the  periodic  temperacure  fluctuations  of  rocacing 
?HF  are  examined.  Finally,  for  different  flow  races  of  the  combustible 
mixture  througn  the  burner,  the  flame  temperature  corresponding  co  cr.e  onset 
of  cell  formation  are  determined. 

EXPERIMENTAL  FROCZDIHES 

The  ?HF  of  rich  butane-air  mixtures  are  scaoiliced  on  the  3uasen  burner 
shown  in  Fig.  1.  The  burner  i3  a  brass  tube  -with  1.2T  and  9.96  cm  outside  and 
inside  diameters  and  is  1  a  long  such  that  fully-developed  parabolic  velocity 
orociles  occur  at  the  burner  exit  olane.  .An  outer  concentric  cube  with  2.J  at 


inside  diameter  is  provided  for  the  control  of  velocicv/ composition  of  the 
surrounding  snosphera.  Four  equally-spaced  layers  of  snail  re sn- screens  are 
placed  within  the  annular  region  between  the  concentric  tubes  to  oocain 
unifora  velocity  profile,  as  shown  in  rig.  2a.  The  fuel,  oxygen,  nitrogen  and 
CF^r  are  separately  aecered  and  may  be  incroaucad  in  the  inner  or  oucar 
tubes.  The  gases  are  thoroughly  nixed  in  a  nixing  chamber  before  introduction 
into  the  burner  tube.  Precise  adjustment  of  the  fuel  and  CF-3r  flow  races  are 
nandacory,  because  small  changes  in  concentration  of  these  species  car  altar 
the  ?HF  substantially  (3idingar  and  Sohrab,  1986) .  Soap-bubble  necers  are 
therefore  used  to  accurately  measure  the  flow  races  of  the  butane,  CF-;3r  and 

A  stationary  6— sided  ?HF  is  chosen  as  a  standard  reference  flame  for  the 
subsequent  comparisons.  For  Che  reference  flame,  the  fuel  concentration  (* 
by  volume)  is  4.3  and  the  mean  velocity  at  the  burner  exit  plane  V  is  72.65 

ea/s.  The  light-blue  flame  with  a  closed  tip  has  the  height  of  about  13  mm, 

and  the  dark-zone  between  the  burner  ria  and  the  flame  base  is  0.3  ma.  In 
order  to  prevent  the  flame  rotation  about  the  axis  of  symmecry,  during  the 
temperature  measurements ,  a  special  stabilizer  ring  shown  in  Pig.  2b  was 
designed.  .The  stabilizer  ring  is  mounted  at  the  top,  around  the  oucar  burner 
tube  such  that  it  can  be  rotated  about  the  burner  axis.  Six  platinum  wires, 

0.05  am  diameter,  are  circumferentially  positioned  such  chac  they  are  60 

degrees  apart  and  extend  radiaily  inward.  The  'wires  terminate  3.5  nm  away 
from  the  inner-surface  of  the  burner  cube,  close  to  the  cresc  zones  of  the 
reference  ?HF.  Rotation  of  the  stabilizer  ring  results  in  rigid-body  rotation 
of  the  encira  ?HF.  It  is  established  that  the  stabilizer  ring  has  no  visible 
effects  on  the  reference  ?HF.  This  ring  helps  to  prevent  the  unwanted 
rotation  of  the  flames  during  the  temperature  measurements  on  stationary 


?KF.  A  direct  pnotograpn  of  tr.e  standard  ?HF  :a  the  ourner  wi : staoil.cer 
ring  is  shown  in  Fig.  3. 

Teraperacura  measurements  ara  made  by  using  a  silica-coated ,  ?t-?t-I0Z  Rh 
thermocouple  vich  0.05  an  wire  diameter,  which  can  be  traversed  in  boch  radial 
and  axial  directions.  The  response  ciae  of  the  thermocouple  junction  is 


estimated 

to  be 

about  2  ms.  3ecause  of  the 

local  perturba 

tions 

induced 

by  t  he 

intrusion 

of  the 

thermocouple,  the  precise 

positioning  of 

the  j 

unction 

is  very 

difficult. 

For 

example,  if  the  junction  is 

positioned  too 

far 

upstream 

of  the 

flame  sheet  the  flame  sheet  aay  corrugate  and  anchor  at  the  junction  . 

However,  .since  the  luminous  flame  sheet  is  clearly  identif iable ,  the 
approximate  positioning  of  the  junction  in  the  vicinity  of  the  reaction  cone 
can  be  achieved.  It  is  found  necessary  to  employ  two  judgemental  criteria  for 
the  positioning  of  the  thermocouple  junction. 

The  following  criteria  are  used  for  the  positioning  of  Che  thermocouple 
junction  in  che  trough  and  crest  tones  of  stationary  ?HF.  The  ?HF  is  anchored 
on  che  burner  and  held  stationary  by  the  stabilizer  ring.  The  ring  is  then 
rotated  and  the  flame  is  positioned  such  that  the  center  of  trough  zone 
radially  points  toward  che  thermocouple  junction.  The  junction  is  then 
radially  moved  inward  in  such  a  way  that  che  junction  just  touches  the 
iownstream  edge  of  the  luminous  cone  of  the  flame  sheec.  The  cri terra  cor 
oosicioning  of  che  junction  in  the  crest  is  different.  Here,  the  junction  is 
positioned  at  the  outer  edge  of  the  dark,  crest  cone  and  is  symmetrically 
Located  between  two  adjacent  calls.  The  adjustment  of  the  position  has  to  be 
lone  as  oreciseiy  as  possible,  because  the  smallest  deviations  may  cause 
rather  large  temperature  variations.  Falsification  of  the  measurements  can 
occur  bv  attachment  of  the  flame  sheets  of  the  adjacent  ceils  to  the  junction 


located  in  the  middle 


The  procedure  described  above  was  developed  such  :na:  300c  repeatability 
of  che  measurements  could  be  achieved.  The  temperatures  of  che  crougn,  T. , 
and  cresn  (ridge),  T,.,  2ones  are  chen  measured  with  che  junction  locaced  5  era 
above  che  burner  rim.  Variations  of  T.  and  T„  are  determined  as  che 

w  «. 

voluraecric  concencracions  of  che  inhibicor,  2 _  and  nitrogen,  p  are 

C7  3r  5  X 

gradually  increased  in  che  reference  PHF  until  complete  disappearance  of  che 
cellular  structure. 

For  conducting  camperature  measurements  in  rocacing  ?Hr,  che 
thermocouple,  an  amplifier,  and  an  oscilloscope  are  arranged  according  co  che 
schematic  drawing  shown  in  Fig.  4.  The  stabilicer  ring  is  rsraovisd  for  tests 
on  rotating  PHF.  Tests  on  rotating  PHF,  consider  butane  concentration  of  4.35 
in  air  at  velocity  7  *  72.65  ca/s  suggested  by  che  earlier  observations 
(Sohrab  and  Law,  1985).  The  introduction  of  nitrogen  into  che  outer  jackec 
Cube  removes  che  surrounding  diffusion-flame  mantle  and  results  into  rapid 
rotation  of  che  PHF.  The  thermocouple  junction  is  located  at  a  constant  axial 
position  6  mm  above  che  burner  rim.  If  the  junction  is  coo  close  c£>  che 
burner  rim  it  will  prevent  the  rotation  of  che  PHF.  Temperature  of  che 
upstream  and  downstream  regions  of  che  rocacing  flame  are  cnen  measured  by 
varying  the  radial  position  of  the  junction  from  2.4  co  4.->  mm  relative  co  che 
axis  of  symmetry. 

In  order  co  decect  che  number  of  cells  of  che  rapidly  rocacing 
polyhedrals,  cwo  thermocouples  are  employed  and  connected  in  series,  as  shown 
in  Fig.  3.  The  variation  of  the  phase  angles  between  the  two  thermocouples  is 
chen  determined  from  che  output  signal  when  one  junction  is  fixed  and  the 
ocher,  at  a  fixed  radial  and  axial  positioa,  is  circumferentially  rotated. 

When  che  signal  of  boch  junctions  are  in  phase,  say  both  locaced  ac  che  crest 
or  trough  position,  che  angle  beeveen  :he  junccion  will  correspond  co  che  cell 


angle.  The  number  of  sides  of  :ne  rocatir.g  ?HF  is  then  data  —  ined  froc  the 
’Ar/v-rlet!«a  of  ceil  angles.  Tor  example,  for  a  o-siced  stationary  ?Hr,  the 

ceil  angle  will  be  60  degrees. 


STATIONARY  POLYHEDRAL  FLAMES 


As  the  additives  C7j3r  and  N^  are  gradually  introduced  into  the  reference 

?HF,  the  cellular  flame  structure  is  found  to  diminish  and  eventually  a 

conical-flame  with  smooth  surface  is  achieved.  The  minimum  values  of  2 _ 

IE,  3r 

and  1  corresponding  to  smooch  flame  situation  are  rasnect ivelv  0.36  and  5* 

■A 

in  close  agreement  with  the  previous  observations  (3idinger  and  Sohrab, 

1986).  Also,  as  CF^Sr  (N^)  is  added,  the  flame  height  increases  from  18  tan  to 

32  mm  (30  am)  at  the  smooch  flame  condition.  Direct  photographs  of  the 

reference  PHF  before  and  after  addition  of  0.363  CF-j3r  are  shown  in  Fig.  6. 

It  is  noted  that  the  difference  between  influences  of  CE33r  and  No  cannoc  be 

accounted  for  on  the  basis  of  thermal  effects  alone.  For  the  typical 

tamoerature  of  1670  K,  the  molar  soecific  heats,  c  ,  of  CF~ 3r  and  N-  are 

?i 

respectively  25.17  and  8.43  cal/mole-'<.  When  account  is  made  of  the 
respective  noia  fractions,  x,  ,  of  these  additives  (3.36  and  53),  the  ratio 
between  x  c  is  about  3.27.  Therefore,  similar  to  the  flame  inhibition 


studies  (Sohrab,  1981)  it  is  concluded  that  CF^r  effects  occur  through 
chemical-kinetic  interference  with  the  flame  structure. 


The  variations  of  T_  and  T_  as  a  function  of  _  and  a,  ara  shown  in 

CE^3r  TI, 

Figs.  7  and  3.  Examination  of  the  data  shows  chat  both  T=  and  Tr  decrease 

with  .1  ,  while  they  first  increase  and  then  decrease  with  .  It  is 

^2  w.--3r 

interesting  that  the  critical  value  of  „  associated  with  the  maximum 

C?33r 


tamperacures  T_  and  T«  in  Fig.  7,  about 


closely  agrees  with  the  critical 


n  i 

W  •  i  , 


inhibitor  -sola  fraction  corresponding  to  me  maximum  inhibition  efficiency  of 
CF^Br  determined  earlier  (Sohrab,  1931).  In  the  previous  study,  inhibition  of 
diffusion  fiames  burning  in  the  stagnation-point  flow  above  liquid  pool  of 
hepcane  was  investigated  (Sohrab,  1981).  It  was  found  chat  the  maximum  oxygen 
mass- fraction  Tq  at  extinction  occurs  when  the  mole  fraction  of  inhibitor  in 
C?3 3r-0-> -N.,  mixture  is  x  ■  0.1  as  shown  in  rig.  9.  Therefore,  the  data  in 
Figs.  7  further  confirms  that  the  inhibition  efficiency  of  CF-3r  is  reiaced  to 
the  temperacura— sensitivities  of  the  kinetics  of  inhibition  as  'was  discussed 
previously  (Sohrab,  1981).  The  comparisons  of  the  data  in  Figs.  7  and  3  also 
show  that  temperatures  T.  and  Tr  are  higher  in  the  flames  inhibited  by  CF ^ 3r 
which  is  a  manifestation  of  the  chemical  rather  than  thermal  influence  of  this 
additive  as  compared  to  nitrogen. 

In  a  series  of  tests,  Che  temperature  of  the  flames  corresponding  to  the 
onset  of  cellular  structure,  TQ ,  is  determined  at  different  gas  velocities 
V.  Here,  onset  of  cellular  instability  is  defined  as  the  first  appearance  of 
a  dark  region  over  a  narrow  vertical  strip  on  the  smooch  conical  flames. 
Starting  wich  smooch  rich  flames,  n_  is  steadily  reduced  'until  the  first 
appearance  of  Che  dark  zone  at  which  point  the  temperature  of  the  luminous 
zone  of  the  flame  is  measured.  The  data,  presented  in  Table  1,  show  that  T3 
is  nearly  constant,  i  10  °C,  for  various  velocities  7.  Because  of  the 
uncertainties  in  the  judgement  concerning  the  initial  appearance  of 
instability,  the  data  are  qualitative  in  nature.  However,  the  results  appear 
to  support  the  fact  chat  chemical-kinetic  mechanisms  may  be  effective  in 
sustenance  or  amplifications  of  the  flame-front  perturbations.  Such  behavior 
is  somewhat  similar  to  the  well  known  constancy  of  the  flame  temperatures, 
about  1120  °C,  at  che  lean  flammability  limits  of  lower  hydrocarbons,  such  as 
methane,  orooane  and  butane.  This  latter  phenomenon  suggests  the  possible 


role  of  chemical-kinetic  mechanisms  in  -he  description  of  the  lean 
flammability  limits  of  the  lighter  hydrocarbons.  Such  chemical-kine tic 
effects  are  further  discussed  in  the  following  section. 

PHENOMENOLOGICAL  ARGUMENTS 

The  importance  of  molecular  hydrogen  in  the  formation  of  cellular  flames 
was  emphasized  in  the  early  reports  by  3ehrens  (1953)  and  Garside  and  Jackson 
(1953).  It  was  found  that  addition  of  small  quantities  of  hydrogen  will 
enhance  the  tendency  far  ?HF  formation  in  propane-air  flames  (Garside  ana 
Jackson,  1953).  In  particular,  no  ?HF  could  be  obtained  in  methane-air 
mixtures  unless  some  hydrogen  was  also  added  (3ehrens»  1953).  Certain 
phenomenological  arguments  were  presented  by  3ehrens  (1953)  in  support  of  his 
experimental  observations.  It  was  also  noted  that  the  chemical-kinetics  of 
methane  combustion  differs  from  that  of  the  higher  hydrocarbons  in  chat  CH^ 
radical  instead  of  the  atomic-hydrogen  plays  an  important  role  in  methane 
combustion.  In  this  pioneering  work  (3ehrens,  1953)  the  author  concluded  tha 
the  molecular  hydrogen,  and  Che  resulting  atomic-hydrogens  as  reaction 
intermediaries,  play  a  key  role  in  the  formation  of  cellular  flames. 

In  the  present  study,  the  significant  role  of  acomic-hydrogan  is  further 
isolated  through  the  introduction  of  the  chemical- inhibitor  CF-;3r.  On  the 
basis  of  the  earlier  ( 3idinger  and  Sohrab,  1986)  and  present  observations,  it 
may  be  concluded  that  chemical-kinetics  effects  have  an  appreciable  influence 
on  cellular  flames.  In  particular,  complete  destruction  of  cellular  flame 
structure  i3  achieved  through  chemical  interaction  of  CF^Sr  with  the  reaction 
zone  of  the  unstable  flame.  Therefore,  the  results  of  the  present  study  are 
better  understood  when  viewed  within  the  framework  of  the  kinetics  of 
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inhibition  by  haioganacad  compounds  reviewed  previously  (Sohrab,  1981).  Here, 
in  Che  same  spiric  as  that  considered  by  Behrens  (1953),  two  phenomenological 
arguments  will  be  presented  which  may  in  part  describe  how  chemical-kinetic 
mechanisms  can  influence  the  flame-front  stability.  The  arguments  are  based 
on  the  kinetics  of  branched-chain  reactions  which  are  relevant  to  hydrocarbon 
combustion  (Lewis  and  von  Elbe,  1967). 

The  central  role  of  atomic  hydrogen  in  the  combustion  of  hydrocarbons  was 
reviewed  previously  (Sohrab,  1981).  It  is  known  that  within  a  chin  radical- 
production  zone,  which  is  embedded  in  the  relatively  wider  reaction  zone,  the 
following  important  reactions  occur: 

a  +  o7  ♦  oh  +  o  (i) 

OH  +  a2  ♦  OH  +■  3  (2) 

OH  +  H2  +  H20  +  a  (3) 

The  race  of  reactions  (2)  and  (3)  are  orders  of  magnitude  higher  than  that  of 

reaction  (1).  Hence,  the  important  chain-branching  reaction  (1)  with 
activation  energy  of  17  kcal/mole  is  rate  controlling. 

The  inhibitor  CFjBr  is  expected  to  first  decompose  according  to  the 
reaction  (Sohrab,  1981), 

CF  3r  *  3  *  C?3  HBr  (4) 

The  subsequent  role  of  CFj  is  not  significant  and  cyclic  removal  of  atomic 
hydrogen  by  HBr  is  expected  to  occur  through  the  regenerative  reaction  cycles 
(Sohrab,  1981). 

3  +  HBr  -  3,  f  3r  (5) 

3r  +  3r  +•  M  *  3r^  +  M  (6) 

and 


Br^  +  H  *  HBr  +  Sr  (7) 

Thus,  the  overall  effect  of  the  inhibitor  is  replacement  of  active  3  acorns  by 
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lass  reactive  bromine  atoms,  thereby  reducing  the  overall  reaction  rate. 

Since  atomic  hydrogen  generation  is  controlled  by  the  reaction  (1)  which 
has  a  large  activation  energy,  the  local  peaks  in  temperature,  say  in  the 
trough  cones  of  PHF,  are  expected  to  result  in  H-atom  concentration  peaks. 
Also,  high  H  concentration  will  enhance  the  consumption  of  C>2  by  the  reactions 
(l)-(3).  This  will  result  in  larger  concentration  gradients  of  0-,  which  will 
in  turn  enhance  the  diffusion  of  O2  towards  the  trough  zone.  Thus, 
stratification  of  the  mixture  could  be  augmented  through,  such  chemical-kinetic 
mechanisms.  As  a  result,  the  temperature  perturbations  which  are  originated 
by  preferential-diffusion  are  amplified  thereby  increasing  the  propensity  for 
cell  formation. 

We  note  that  because  of  the  nature  of  branched-chaia  reactions,  only 
small  deviation  of  La  from  unity  are  needed  in  order  to  produce  the  original 
temperature  corrugations  whose  subsequent  amplification  leads  to  large 
observable  corrugations  of  the  front.  Also,  in  accordance  with  the 
observations,  the  amplification  of  corrugations  is  expected  to  be  reduced  in 
presence  of  since  the  inhibitor  removes  the  atomic-hydrogen  radicals. 

A  schematic  drawing  of  the  expected  fluctuations  of  T,  a,  0-,  and  fuel  F,  for  a 

6-sided  PHF  is  shown  in  Fig.  10.  Such  fluctuations  of  0-,  and  F  are  in 

qualitative  agreement  with  the.  observations  of  Jost,  et  al.  (1953).  Also,  the 
variation  of  T  is  in  accordance  with  the  observations  made  herein  as  well  as 
the  theoretical  predictions  ( Sivashinsky ,  1976).  The  schematic  concentration 
■profiles  in  Fig.  10  are  also  supported  by  the  previous  observations  of  soot 

formation  in  the  crest  zones  of  certain  PHF  (Smith  and  Pickering,  1929).  In 

the  crests,  CO  and  Oj  concentrations  are  larger  (Jost,  et  al.,  1953)  while  H, 
and  therefore  OH,  concentrations  as  well  as  the  temperature  are  lower  as 
compared  with  the  troughs.  Therefore,  formation  of  soot  in  the  crests  is 


anticipated  on  the  basis  of  high  concentration  of  CO,  possible  catalytic 
effects  of  0-,  (Lin  and  Sohrab,  1986)  as  veil  as  the  veil  known  role  of  OH  in 
the  combustion  of  soot  pracursers. 

The  second  phenomenological  argument  concerning  the  role  of  chemical- 
kinetics  in  flame  instability  is  based  on  the  concept  of  the  critical  cross¬ 
over  temperature  (Sohrab,  1981).  3ecause  of  their  different  activation 
energies,  the  rates  of  bimolecular  and  teraolecular  reactions  may  cross  each 
ocher  when  ploctad  as  functions  of  1/T.  The  cross-over  point  becveen  the  rat a 
of  an  important  biomolacular  reaction,  such  as  reaction  (1),  and  an  important 
teraolecular  reaction,  such  as  reaction 

H  +  CL  +  M  ♦  HO,  +  M  ,  (3) 

can  be  defined  as  the  critical  cross-over  temperature  Tc  as  shown  in  rig.  11. 
Vhen  cemneracure  perturbations  about  T  occur,  for  T  >  T  and  T  <  T  the 
reaction  rate  is  expected  to  respectively  accelerate  or  decelerate.  This  will 
result  in  corrugation  of  the  flame  front  as  shown  in  Fig.  II.  Again,  because 
of  the  nature  of  branched-chain  reactions,  only  small  perturbations  about  Tc 
are  necessary  to  produce  substantial  change  in  the  reaction  rate  leading  to 
formation  of  cellular  structure.  The  small  temperature  fluctuations  about  T, 
could  be  initiated  by  the  pref erential-dif fusion  mechanism.  Such  a 
phenomenological  description  is  supported  by  the  constancy  of  temperature  T0 
corresponding  to  the  onset  of  cell  formation  mentioned  in  the  previous 
section. 

ROTATING  POLYHEDRAL  FLAMES 


In  this  section,  the  periodic  variation  of  the  temperature  field  in  the 
vicinity  of  rapidly  rotating  ?HE  is  discussed.  The  fuel  concentration  of  the 
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reierence  ?HF  is  now  changed  to  Cl_  ■  i.95,  with  V  a  72.67  cm/  s  and  che 
stabilizer  ring  is  removed.  The  resulting  stationary  6-sided  PHF  will  start  a 
rapid  counter-clockwise  rotation  when  small  flow,  72  cn^/s,  of  nitrogen  is 
introduced  in  the  outer  annular  tube  as  discussed  before.  A  direct  photograph 
of  Che  rotating  PHF  is  shown  in  Fig.  12.  The  luminous  zone  of  the  rapidly 
rocating  PHF  is  fuzzy  and  thicker  while  the  flame  height  is  smaller  compared 
to  the  stationary  PHF.  Also,  che  rotation  of  PHF  is  found  to  predominantly 
occur  in  che  counter-clockwise  direction  in  agreement  with  che  previous 
observations  (Sohrab  and  Law,  1985). 

The  number  of  sides  and  che  rotation  speed  of  che  PHF  are  measured  using 
che  two— Chermcouple  connections,  Fig.  5,  discussed  earlier.  It  is  found  that 
in-phase  signal  of  che  two  thermocouples  are  obcained  when  che  second  junction 
Is  located  at  90°,  130°  and  270°  angles  relative  to  che  first  junction. 
Therefore,  the  original  stationary  6-sided  PHF  has  changed  to  a  4-sided 
structure  after  rotation  Is  Induced  by  nitrogen  flow  in  the  outer  annulus. 

This  result  Is  in  agreement  with  che  previous  observation  of  rapidly  rocating 
PHF  using  high  speed  cinematography  (Sohrab  and  Law,  1985).  Also,  from  the 
period  of  temperature  fluctuations,  32  ms,  che  rotacion  speed,  w  *  7.8 
rotations  per  second  (rps),  is  calculated  In  qualitative  agreement  with  the 
previous  observation  (Sohrab  and  Law,  1985). 

The  circumferential  temperature  fluctuations  at  different  radial 
positions  are  obcained  from  the  tracing  of  the  direct  phocographs  of  che 
oscilloscope  screen  and  are  shown  in  Fig.  13a-t.  The  variations  of  periodic 
temperature  profiles  at  different  radial  positions  are  qualitatively 
understood  on  the  basis  of  the  residence  time  of  the  thermocouple  junction  in 
che  downstream  (hoc)  or  upstream  (cold)  sides  of  the  rotating  PHF.  In 
general,  che  trough  zones  are  circumferentially  wider  than  che  crest  zones,  as 


reflected  by  the  longer  time  Intervals  of  high  versus  low  temperatures  in  rig. 
13.  Also,  the  magnitude  of  the  temperature  fluctuations  are  a  function  of  the 
relacive  instantaneous  position  of  the  junction  with  respect  to  the  corrugated 
flame-front.  We  note  that  the  temperature  gradients  downstream  of  the  flame 
are  much  smaller  than  those  in  the  upstream  region,  as  is  to  be  expected.  The 
results  in  rig.  13  also  show  that  temperature  fluctuations  occur  about  a  mean 
temperature  plateau  of  approximately  1408  °C.  Hence,  rotation  of  different 
cells  of  the  PHF  occur  in  preheated  gas.  Is  is  known  that  the  rotation  of  the 
flame,  reaction  zone,  is  not  accompanied  with  any  appreciable  circumferential 
rotation  of  the  combustible  gas  (Smith  and  Pickering,  1929;  Sohrab  and  Law, 
1985).  ' 

Rotating  PHF  was  also  examined  in  oxygen-enriched  air.  For  the  fuel  and 

oxygen  concentrations  of  ft„  ■  7.76  and  ft  ■  26. a  stationary  8-sided  PHF  is 

J?  o 

obtained  with  well  defined  cells  and  a  scream  of  soot  is  observed  to  emanate 

from  the  flame  tip.  However,  this  flame  was  found  to  be  only  marginally 

stable,  and  under  small  laboratory  fluctuations,  it  could  spontaneously  start 

a  rapid  counter-clockwise  rotation.  We  note  chat  3uch  rapidly  rotating  PHF  in 

butane/ air  mixtures  could  only  occur  in  the  presence  of  nitrogen  flow  in  che 

outer  annulus  (Sohrab  and  Law,  1985).  Therefore,  the  present  observation 

shows  that  rapidly  rotating  flames  may  occur  in  oxygen-rich  air  in  absence  of 

inert  surrounding  atmosphere.  Here,  the  period  of  temperature  fluctuation  is 

20  ms,  and  the  PHF  is  found  to  be  4-sided,  resulting  in  w  *  12.5  rps . 

According  to  some  preliminary  observations,  Che  rotation  speed  of  the  flame 

was  found  to  increase  with  ft  . 

o 


CONCLUDING  REMARKS 


Through  comparisons  of  che  minimum  mole  fractions  of  che  additives  CE-jBr 
and  N^  required  for  destruction  of  cellular  structure  of  butane/air  ?HF,  the 
significance  of  chemical  versus  thermal  affects  in  the  formation/ amplification 
of  cellular  flame  structures  was  investigated.  The  small  changes  in  che  mean 
specific  heats  and  velocities  induced  by  the  additives,  on  che  ocher  hand, 
revealed  chat  che  observed  cell  destruction  does  not  appear  to  have  thermal  or 
hydrodynamic  origins.  It  was  therefore  concluded  chat  the  influence  of  CF^Br 
primarily  occurs  through  chemical-kinetic  modification  of  che  structure  of 
?HF.  It  is  emphasized  chat  che  kinetic  aechanicsms  discussed  herein  are 
complimentary  to  che  well-known  thermo-diffusive,  Le  ?•  1 ,  mechanisms  of  cell 
formation.  Indeed,  che  observations  are  consistent  with  both  mechanisms 
working  in  conjunction  with  each  ocher. 

Since  pyrolysis  of  fuel  molecules  occurs  prior  to  their  consumption  in 

the  reaction  zone,  more  Chan  one  Lewis  number,  Le^  *  a/D^ ,  each  corresponding 

to  che  diffusion  coefficient  D*  of  specie  i,  could  effect  the  flame  front. 

Here,  a  refers  to  che  thermal  diffusivicy  of  the  mixture.  Therefore,  in 

general,  thermo-diffusive  description  of  flame-front  instability  could  be 

based  on  che  deviation  of  Le  of  reaction  intermediaries,  such  as  radical 

+  * 

species,  from  unity.  In  fact,  recent  theoretical  scudy  of  Pelaez  and  Linan 
(1985)  has  shown  that  thermo-diffusive  stability  of  flames  are  influenced  by 
the  Lewis  number  of  intermediate  species.  •  Sdch  effects  are  particularly 
relevant  to  atomic  hydrogen  because  of  its  high  mobility.  However,  che  role 
of  chemistry  is  separate  and  fundamentally  different  from  that  of  the 
transport,  Le  ^  I ,  effects.  The  former  may  manifest  itself  through 
augmentation  of  preferential-diffusion  effects  by  chemical-kinetic  mechanisms 
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as  discussed  earlier. 


Intrusive  temperature  measurements  in  meeting  ?HF  proved  to  be 
difficult.  However,  the  technique  of  employing  two  thermocouple  junctions 
connected  in  series  facilitated  the  determination  of  the  number  of  cells  as 
well  as  the  rotation  speed  of  PHF.  So  far,  the  reason  for  predominant 
rotation  of  ?HF  in  counter-clockwise  direction  remains  a  mystery.  The 
preliminary  tests  on  the  influence  of  applied  magnetic  field  on  rotating  ?HF 
showed  no  observable  effects.  Thus,  influences  associated  with  the  asymmetry 
in  the  Earth's  magnetic  field  were  ruled  out.  Exploration  of  the  possible 
effects  of  coriolis  acceleration  on  this  interesting  phenomena  through 
experiments  on  rotating  ?HF  performed  In  Che  Souchern  hemisphere,  may  prove 
useful. 
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Abstract— In  this  paper  we  consider  the  influence  of  centrifugal  and  Coriolis  accelerations  on  the  form 
and  stability  of  a  flame  propagating  in  a  rotating  tube.  We  detect  a  discrete  sequence  of  angular 
velocities,  near  which  the  flame  propagation  velocity  undergoes  an  unlimited  amplification.  We  show 
that  sufficiently  rapid  rotation  may  suppress  cellular  instability  of  the  flame  front.  We  construct  a  self¬ 
similar  solution  that  describes  a  flame  propagating  in  a  rotating  gas  in  free  space.  Our  predictions  of  the 
flame  shape  and  of  the  suppression  of  cellular  instability  in  the  presence  of  rotation  have  been  verified 
in  a  simple  experiment  vve  nave  performed. 


INTRODUCTION 

It  is  well  known  that  variations  in  the  density  of  a  burning  gas  make  the  latter 
sensitive  to  external  acceleration  (see.  e.g.,  Rakib  and  Sivashinskv.  1 987 i.  In  the 
present  paper  vve  shall  investigate  the  influence  of  centrifugal  and  Coriolis 
accelerations  on  the  form  and  stability  of  a  premixed  flame  propagating  in  a 
rotating  cylindrical  tube. 

In  order  to  capture  the  nature  of  the  rotation-induced  effects,  we  begin  with  the 
simplest  hydrodynamic  model,  considering  the  flame  to  be  a  geometric  surface 
moving  at  constant  velocity  Uh  relative  to  the  burnt  gas.  Transport  and  chemical 
effects  are  ignored,  but  the  change  in  gas  density  across  the  flame  front  is  taken 
into  consideration  (Landau.  1944). 

In  a  coordinate  frame  rotating  together  with  the  tube  about  the  ti-axis.  at  angular 
velocity  a>.  the  Euler  equations  are: 

dv,  dv,  if  dv.  dvr  if  _  1  dp  :  ^  ^ 

dt  dr  r  Orp  dz  r  a  dr 

dv,r  dv,,  v.,  dv,t  dv,.  i\v,T  1  dp 

— —  +  vr— — i - - — rV-~Z — 1 - - - r - 2u>vr ,  ,1.2’ 

dt  dr  r  dep  '  dz  r  pr  dep 

dv.  dv.  v,r  dv.  '  dv.  _  1  dp 

dt  dr  r  dep  '  dz  p  dz 

dvr  1  dv.  dv.  vr 

T-  +  -“  +  ^+-  =  °.  '1.4' 

dr  r  dep  dz  r 


Here  [r,  rp.  z)  are  nondimensional  cylindrical  coordinates,  with  r,  c  measured  in 
units  of  the  radius  R  of  the  tube  (see  Figure  1);  t  is  the  nondimensional  time, 
measured  in  units  of  R/U,,\  to  is  the  nondimensional  aneular  velocity,  measured  in 
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units  of  Ub.R\  o  is  the  nonciimensional  density  measured  in  units  of  the  density  nh 
of  the  combustion  products;  p  is  the  nondimensionai  pressure,  measured  in  units 
of  ohUh:  and  (v„  v r  v.)  =  v  is  the  nondimensionai  velocity  of  the  gas,  measured  in 
units  of  Uh. 


flame 


FIGURE  1  Schematic  representation  of  a  flame  propagating  in  a  rotating  channel. 


The  following  conditions  must  hold  at  the  flame  front  [z-  <I>(r,  <p,  /)]: 

i)  Continuity  of  mass  flow’ 

[p(v-n-  D)}  =  0, 

ii )  Continuity  of  momentum  flow 

|/)v(vn-£))+/tnj  =  0.  ;l-6) 


iii)  Constant  velocity  of  the  flame  front  relative  to  the  gas 

,oiv-n  —  D}  =  l.  1 1.7) 


where 

*„  l| 

'  '  !  1.S1 

V1  +,VcJ>  r  ’  vl+(V<D,r 

Hydrodynamic  quantities  corresponding  to  the  burnt  gas  region  [t;^<t>(r,  cp.  r)]  will 
be  assigned  the  index  /  - );  those  corresponding  to  the  fresh  gas  region 
■  C  <  <I>(  r~tp,  />]  the  index  t  -  ).  With  this  convention. 

,o  +  ~  1 ,  p-  =(1  ~  yf\  !l-9) 


SELECTION  OF  SCALINGS 

In  the  absence  of  rotation  (a>  =  0).  one  of  the  possible  states  of  the  system  is  a 
plane  i  p\  flame: 

=  0.  vi  =  (0.  0,  1 ),  vi  =  1 0.  U.  1  -  yj,  pi  =  q'.  —  1 .  pi  =  qi  =  1  -  y.  2.1 ; 

If  the  flame  is  propagating  in  a  rotating  tube,  the  plane  flame  can  no  longer  be  an 
equilibrium  state  of  the  system.  Rotation-induced  acceleration  produces  a  radial 
flow  which  curves  the  flame  front.  However,  if  the  thermal  expansion  of  the  gas  is 
weak  t  1 )  the  curving  will  be  slight.  Thus,  one  expects  the  equation  for  the  flame 
front  dynamics  to  be  weakly  nonlinear.  The  structure  of  this  equation  is  largely 
determined  by  condition  11.7),  which  assumes  a  simpler  form  when  y  <  1 : 

<p  +1 ,  V<t>  ::  =  v‘.  -  l.  2.2' 

Condition  (1.15)  suggests  that  when  o>=  0(1)  the  perturbation  of  the  reduced 
pressure  q  due  to  rotation  is  of  the  order  of  magnitude  of  y.  It  is  natural  to  assume 
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that  the  perturbation  of  the  velocity  v  due  to  rotation  will  be  of  the  same  order  of 
magnitude.  According  to  this  estimate  and  {2.1;,  we  introduce  scaled  qualities  l,’, 

v\  —  yV~, ,  = 

v]  =  l  +  yV+.,  v~.  =  l-y+yV~.  , 

q*  =  1  -  O",  c/==l  +y+yQ~-  2.3) 

Since  r~  1  in  terms  of  the  selected  nondimensional  units,  the  relationships  2.2) 
and  (2.3)  suggest  the  following  scalings  for  the  time  /and  the  front  amplitude  (t>: 

<$>  =  JyF,  t=x/]y.  (2.4) 


In  terms  of  the  scaled  variables  (2.3),  (2.4),  we  obtain  from  (1.1)—'  1.4)  the 
following  equations  for  the  principal  term  of  the  asymptotic  expansion  i  y  <  1 ): 


a  i/-;  dQz 
dz  dr 


2coV;r. 


dVl  ,  1  dQ'~ 
dz  r  dcp 


-2 cjV:  . 


dV :  dQz 
dz  dz 


=  0, 


(2.6) 


(2.7) 


dr  r  dcp  dz  r 


/  2.8 ) 


Substituting  *  =  =  v'yF  in  condition  ( 1.1 1 W  1.13)  and  (1.7)  and  taking  the  zeroth 

approximation  with  respect  to  y,  we  obtain  the  following  conditions  at  z=  0: 


v;  =  vr ,  v:  =  •; ,  v:  =  v: 


2.9' 


Q~  -  Q  =  G<> i 


-<u'r , 


2.10) 


I/:  .  .2.11; 

Impermeability  conditions  are  imposed  a'  the  tube  walls  (r=  1 ): 

K  =  v:=0  at  r=l.  (2.12) 

We  shall  also  assume  that,  far  ahead  of  the  flame  front,  the  gas  rotates  as  a  rigid 
body.  Thus,  in  a  rotating  coordinate  system,  we  have  the  following  conditions  as 
z  -  -  00 : 

k;~o,  v;-o.  v:-o  as  ?2.i3) 


FLAME  PROPAGATION  IN  ROTATING  GAS 


derivation  of  the  flame-front  equation 

In  the  limit  considered  above,  the  flame  from  configuration  [:=  F' r.  p.  rjj  appears 
only  in  condition  (2.11).  Hence  the  gas  flow  is  completely  determined  by  Eqs. 
(2.5)-(2.8)  and  conditions  (2.9),  (2.10).  (2.12)  and  (2.13)f.  Since  the  problem 
posesses  circular  symmetry,  the  flow  is  independent  of  the  angular  coordinate  cp. 
Hence,  we  seek  a  solution  of  the  linear  system  (2.5 }— (2.8 )  in  the  form 

» 

V.=  1  Lie,  k„;J,){k,,ri, 


K=  I  k„Vr'\z.  kn)Uk„r\, , 

r/-  1 

TO 

K-=  1  k„  Vr (i.  k,,)J„\k„r), 


0=  X  k„}llk„r), 

'f-  1 


where,  by  condition  (2.12), 


X(At„)  =  0.  3.2 

We  also  express  the  right-hand  side  of  (2.10)  as  a  Fourier-Bessel  expansion: 


Gir)=  K  G[k„)J,&k„r), 

n-  I 


J.J  I 


where 


J  A..,, 


rG{r)Jl){k„r)dr= 


.co 


Ji)(  k.,,k„ 


3.41 


Inserting  (3.1)  into  the  system  i2.5>-;2.8)  and  the  conditions  2.9).  .2.101,  we 
obtain: 


dV: 

dz 

dV 


^  Q1  -  -co  , 


u<. 


*  =  -2cuI/7  , 


coai 


3.5b) 


tThe  independence  of  the  hydrodynamic  field  on  the  name  front  configuration  completely  excludes 
the  effect  of  sponateous  name  instability  due  to  thermal  expansion  (Landau.  1 944),  which  turns  out  here 
to  be  negligibly  small  compared  with  the  effects  of  rotation. 
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FLAME  PROPAGATION  IN  ROTATING  GAS 


Hence 


,  ^  r~  Cnk„)k„ 

hT  0;  =  -  _  773 — — y \Uk„r<. 

2  ki  >iu,'.(k.,-4aj  1 


3.13 


From  (2.1 1 )  and  (3.13)  we  obtain  a  single  equation  for  the  flame  front  dynamics 


Fr~- —  iVf '■=  win. 


3.1  4, 


where 


W(r)  =  I  -^4^. 

4u^<it •  J'A  k „ /[  k rl  4  uj  f 

and  use  have  been  made  of  (3.41  and  (3. 13). 

Knowing  W(r<,  we  look  for  solutions  in  a  form  of  a  progressive  wave 

F"  =  -  IT  —  y/(n 


3.15 


5. 1  6  1 


where  V  and  i//  are  the  equilibrium  velocity  and  shape  of  the  flame  front, 
respectively.  Substituting  into  \ 3. 14)  one  obtains 


~  yj]=  W(r)+ V.  3. re 

In  order  for  the  solution  xjir)  to  exist  over  the  entire  tube  cross-section  .0  <  r<  1. 
the  velocity  V  must  satisfy  the  inequality 

V>  -  \V(r).  3.1 8^ 

We  assume  that  the  actually  realized  situation  corresponds  to  minimum  flame 
propagation  velocity,  i.e.. 


V=  -  min  \V(r  1,  \0 <r<  1;. 

The  need  tor  this  selection  principle  arises  here  because  of  the  lack  of  natural 
boundary  condition  for  the  first  order  differential  equation  (3.14).  This  difficulty  is 
easily  eliminated  if  one  introduces  dissipative  effects  diffusion  and  heat 
conduction)  with  the  appropriate  boundary  conditions.  Figures  2  and  3  exhibit  the 
velocity  V  and  the  shape  ip  of  the  flame  front  as  functions  of  the  parameter  4<y:. 

We  note  that  V  does  not  vary  monotonicaily  with  increasing  to.  As  4 uF. 
approaches  k:n,  the  flame  velocity  undergoes  an  unlimited  amplification.  Near  these 
resonance  points  the  asymptotic  approach  proposed  abqve  is  clearly  not 
applicable.  Description  of  the  flame  in  the  vicinity  of  these  points  requires  special 
consideration.  However,  the  non-monotonic  behavior  of  the  flame  with  increasing 
angular  velocity  is  apparently  captured  correctly  by  our  analysis.  It  would  be 
interesting  to  conduct  a  systematic  laboratory  investigation  of  this  phenomenon.  As 
yet,  to  the  best  of  our  knowledge,  the  only  experimental  study  in  this  area  is  the 
recent  work  of  Chen,  Liu  and  Sorhab  (1987),  devoted  to  the  influence  of  rotation 
on  counterflow  premixed  flames. 
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FIGURE  2  Flame  velocity  Kas  a  function  of  the  parameter  4or. 


THE  CASE  OF  SLOW  ROTATION 


If  4ar  <  k,:.  it  may  be  shown  from  (3.15)  that  W(r)  is  a  solution  of  the  equation 

d2W  1  dW  , 

— —  -i - =  a/  , 

dr"  r  dr 


satisfying  the  additional  conditions 


<m 0) 

dr 


=  n 


'  lA/[  /Vr  = 


Hence 
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FIGURE  3  Flame  from  confieuraiion  for  various  values  of  the  parameter  iar  ■  W k:.,  +  k\, 
\1=  1.2.3.  A). 


If  ij>-  <  A:,,  this  relationship  simplifies  to 

Win  —  -  u)'\  r  — j 

*  \  -I 

We  then  deduce  from  1 3.141,  p.i  7  and  (3.18)  that 


4. 


a 


where  V=  ar/S. 

The  stability  of  this  solution  will  now  be  examined. 

For  a  small  perturbation  /=  F-  Eq.  f  3.14)  yields  a  linear  equation 

fr  +  2jVrf  =  0.  f[r.  0 )=Qr\ 


4. 
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Here  jinr;  is  the  initial  perturbation.  It  follows  from  -  4.0  >  that 

fir,  -4.7 

Thus,  small  perturbations  are  stretched  and  the  flame  front  becomes  smooth:  the 
situation  being  similar  to  the  one  found  in  stagnation-point  flow  flames  (Sivashinskv 
er  a/.,  1982),  freely  expanding  flames  (Sivashinskv,  1983)  and  upward  propagating 
flames  (Rakib  and  Sivashinskv,  1987). 

It  should  be  noted  that  when  4 or  <?  k] :  the  effect  of  Coriolis  acceleration  turns 
out  to  be  small  compared  with  the  effects  of  centrifugal  acceleration,  and  the 
corresponding  terms  in  Eqs.  ( 1 . 1  U  J  .2 )  may  be  dropped. 


THE  EFFECT  OF  ROTATION  ON  CELLULAR  INSTABILITY 


Centrifugal  acceleration  creates  a  gradient  of  the  tangential  velocity  which  tends  to 
smooth  out  any  possible  perturbations  of  the  flame  from.  Thus,  one  may  expect 
rotation  to  suppress  diffusional-thermal  i  cellular  i  flame  instability,  frequently 
observed  when  the  deficient  reactant  is  of  high  mobility.  We  shall  demonstrate  this 
for  the  simple  case  of  slow  rotation  considered  in  the  previous  section.  As  a  first 
step,  let  us  write  the  equation  of  the  flame  front  (3.14),  1 4.4 ),  in  terms  of 
dimensional  variables  denoted  by  (*): 


d<t>* 

dr* 


J. 

i 


(A-;v  <t>*r 


i 


- 1 


yQ'/cu; 


where  Q  =  Uhco/R  is  the  dimensional  angular  velocity  as  noted  earlier.  Allowing  for 
diffusional-thermal  effects  implies  that  these  equations  must  be  modified  as  follows 
(see  Sivashinskv,  1983): 


d<&*  1  ,  ; 

— +  -  l/*(V<D*)'  -i-  eD„V-<J>* 
at  2 


4A/i/,/,V4$#  =  |  yOrtflrt 

4 


where  D.h  is  the  thermal  diffusivitv  of  the  gaseous  mixture.  i.h  is  the  thermal 
thickness  of  the  flame  D,h  =  UJ,h/  and  e  =  >  E/2 R" Th,{  1  -  Le:  -  1 .  with  E  the 
activation  energy,  R'  the  universal  gas  constant,  77  the  adiabatic  temperature  of  the 
combustion  products.  Le=  D:/rOmtn  the  Lewis  number  and  Dmol  the  molecular 
diffusivitv  of  the  deficient  reactant. 

Transforming  to  scaled  nondimensional  variables  F.  r.  t 

<D  *  =  ?/,/.  r*=  ZLr/fe.  t*~4l!hi/e:Uh,  5.3' 

we  bring  Eq.  (5.2'  to  the  following  form: 

F.-\vF  :  -  V  :f-*-  V  4F=  Ai'r:  5.4 

where  {•)  denotes  the  spatial  average  over  the  cross-section  and 


k  —  4yc  Uif 


is  the  centrifugal  acceleration  parameter. 
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Figure  4(a,  b,  c)  exhibits  the  results  of  a  numerical  solution  of  the  axiallv- 
symmetric  version  of  Eq.  (5.4)  in  a  tube  0  <  r <  50,  with  boundary  conditions 

£=0,  (V'f);=  0,  at  r=  50,  (5.6) 

corresponding  to  the  thermally  insulated  tube  walls  (Margolis  and  Sivashinskv, 
1984).  At  A  =  0,  Le.,  a  freely  propagating  flame  in  a  nonrotating  channel,  since  e>  0 
the  flame  exhibits  thermo-diffusive  instability  which  manifest  itself  in  the  form  of  a 
wavy  (cellular)  surface  (Figure  4a).  At  A  =  0.001,  cellular  structure  is  superimposed 
on  the  paraboloid  profile  induced  by  rotation  (Figure  4b).  Finally,  when  A  =  0.005, 
centrifugal  acceleration  totally  suppresses  the  cellular  structure  (Figure  4c).  The 
flame  exhibits  a  smooth,  paraboloid  profile,  as  in  the  absence  of  diffusional-thermal 
instability,  i.e.,  when  £<  0. 


ig.4  (a) 
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EXPERIMENTAL  OBSERVATIONS 

The  evolution  of  the  flame  front  shown  in  Figure  4  is  substantiated  bv  our 
preliminary  experiments  performed  on  a  rotating  Bunsen  burner.  In  the 
experiment,  a  cellularly  unstable,  polyhedral  flame  in  rich  butane/air  mixture  with 
volumeric  fuel  concentration  of  6.6  percent  is  stabilized  at  the  exit  plane  of  the 
burner.  The  burner  tube  is  made  of  pyrex  and  is  15  cm  long  with  ir  .de  diameter  of 
2.0  cm,  and  the  mean  axial  velocity  of  the  gas  over  the  cross-section  is  44  cm/s. 
The  tube  is  co-axially  attached  to  a  cylindrical  burner  which  contains  layers  of 
small  mesh  screens  and  honeycomus  and  can  be  rotated  by  a  DC-motor  as 
described  previously  (Chen  et  al.,  1987).  As  the  angular  velocity  Q  is  increased,  the 
cellular  structure  of  the  flame  is  observed  to  grad-ally  diminish.  At  larger  Q. 
Q  —  9  rps,  the  central  position  of  the  flame  cone  is  flattened  until  it  eventually 
becomes  convex  towards  the  unburnt  gas.  At  still  larger  Q,  Q=  12  rps,  the  flame 
moves  into  the  rotating  tube  and  remains  almost  stationary  near  the  middle,  about 
8  cm  below  the  rim,  and  appears  as  an  inverted  smooth  paraboloid.  A  direct 
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Fig.  4(c) 


FIGURE  4  Flame  front  configuration  for  various  values  of  the  parameter  k. 

photograph  of  this  flame  is  shown  ;n  Figure  5.  Qualitatively  similar  observations 
were  made  for  lean  butane/air  mixtures  with  the  exception  that  the  initial  cellular 
structure  was  absent.  The  more  complex  peculiarities  predicted  for  still  larger 
values  of  Q  are  as  yet  not  observed  because  of  limitation  of  the  maximum  angular 
velocity,  Qmax  =  29rps,  of  the  DC-motor.  Further  studies  of  this  interesting 
phenomena  and  experiments  at  larger  Q  will  be  considered  in  the  future. 

FLAME  PROPAGATION  IN  FREE  SPACE 

Hitherto  we  have  been  considering  a  flame  propagating  in  a  rotating  tube.  Under 
these  conditions,  the  flame  front  is  curved.  It  turns  out  that,  if  the  walls  are 
removed,  a  flame  propagating  in  a  rotating  gas  may  remain  planar.  Moreover,  the 
corresponding  gas  flow  is  described  by  a  self-similar  solution  of  the  Euler 
equations,  which  can  be  constructed  with  no  restrictions  imposed  on  the  parameter 
Y- 
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FIGURE  5  Rich  buiane/air  flame  in  a  rotating  tube. 
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Let  us  assume  that  the  flow  in  the  region  of  the  fresh  mixture 
coordinate  frame)  is  unidirectional: 

(in  a  rotating 

u'  =  l-y,  v,=  0,  u"=0, 

•£> 

i 

li 

i 

1 

N> 

-  e 

•  c 

>i 

(6.1) 

In  the  bumt-gas  region,  we  seek  a  solution  of  the  following  form: 

v t  =  2 a[z),  v*  =  -  d[z\r ,  d+,  =  b{z)r. 

*  t  ^  2  2  r"  2  ,  /  \ 

?  2“^“2(l-v)^+cU)- 

(6.2) 

At  the  flame  front,  which  is  assumed  to  be  plane  {z~  O  =  0),  we  have  the  following 
conditions,  derived  from  (1.5),  (1.11),  (1.12)  and  (1.15): 
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II 

►— * 
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II 
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II 
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(6.3) 

Inserting  (6.2)  into  the  Euler  equations  ( 1 . 1 )-( 1 .4 ),  ( 1 . 1 4 ),  we  obtain 

ab'  -  db-  (x)d. 

(6.4) 

2  aa  [a)  +  b  +  2  cub—  , 

1  -y 

(6.5) 

4aa'  +  d  -  0. 

(6.6) 

Note  that  here  the  continuity  equation  ;  natically  satisfied. 

From  (6.4),  (6.6)  and  conditions  '  ,  obtain 

b\z)=2a>a{z)-  a). 

>6.7) 

c{z)  =  ^-2a2{z). 

(6.8) 

(6.10) 
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The  streamlines  corresponding  to  the  flow  (6.2)  are  spirals  winding  around  the 
surface  of  rotation  [see  Eq.  ( 1 .4)]. 


const 


(6.11) 


One  of  these  streamlines  and  the  corresponding  surface  (6.11)  are  illustrated  in 
Figure  6. 


CONCLUDING  REMARKS 

Using  the  weak  thermal  expansion  approximation  we  derived  an  equation 
descnbing  the  dynamics  of  a  flame  propagating  in  a  rotating  tube.  The  velocity  and 
shape  of  the  equilibrium  flame  front  were  calculated.  A  discrete  sequence  of 
angular  velocities  are  predicted  near  which  the  flame  propagation  velocity 
undergoes  a  considerable  amplification.  When  the  angular  velocity  of  rotation  is 
sufficiently  high,  we  predict  that  cellular  structure  due  to  thermo-diffusive  flame 
instability  may  be  suppressed.  The  variation  of  the  flame  shape  and  the  suppression 
of  instabilities  as  a  result  of  the  rotation,  were  experimentally  substantiated.  We 
finally  remark  that,  although  the  analysis  assumes  weak  thermal  expansion  and  in 
real  flames  density  variations  are  ^uite  appreciable,  this  limit  may  give  a 
qualitatively  correct  description  of  some  features  of  real  flames. 
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Abstract 

The  influence  of  centrifugal  accelerations  on  the  form  and  stability  of  the  Bunsen  burner 
stabilized  flame  is  studied.  It  is  shown  that  as  the  angular  velocity  of  rotation  increases,  the 
classical  Bunsen  cone  first  slightly  buckles  at  the  tip.  Buckling  then  becomes  more  prominent  and 
eventually  the  flame  enters  the  burner.  At  sufficiently  rapid  rotation  flame  stabilization  by  the  burner 
rim  becomes  impossible  and  the  flame  flashes  back.  The  theoretical  findings  are  in  favorable 
agreement  with  the  experimental  observations  made  on  a  rotating  Bunsen  burner. 


‘Department  of  Applied  Mathematics,  Tel-Aviv  University,  Ramat  Aviv.  Israel. 


I.  INTRODUCTION 


Flame  propagation  in  situations  of  practical  interest  occurs  as  a  rule,  in  turbulent  flow  fields  and 
involves  complex  flame-flow  interactions.  One  of  the  basic  interactions  of  this  type  is  that  which  is 
imposed  by  rotation  of  the  gas  flow  crossing  the  flame  front. 

Recently,  studying  the  influence  of  centrifugal  and  Coriolis  accelerations  on  the  shape,  stability 
and  extinction  limits  of  premixed  flames  [1]  (cf.  also  [2-5]),  a  number  of  new  observations  were 
made  on  a  rotating  Bunsen  burner.  Among  other  observations,  it  was  found  that  flame  stabilization 
by  the  burner  rim  is  possible  only  at  sufficiently  low  angular  velocities.  Under  rapid  rotation,  the 
flame  flashes  back  inside  the  burner  tube. 

As  we  show  in  this  paper,  the  main  features  of  the  phenomena  are  describable  within  the 
framework  of  a  simple  mathematical  model  recently  proposed  in  the  study  of  premixed  flames  freely 
propagating  in  rotating  tubes  [5]. 


II.  EXPERIMENTAL  OBSERVATIONS 

Premixed  flames  are  stabilized  at  the  rim  of  a  50  cm  long  pyrex  tube  with  11  (13)  mm  inner  (outer) 
diameter.  The  tube  is  vertically  supported  at  20  and  30  cm  from  the  rim  and  at  the  bottom  by 
bearings  and  connected  through  a  pulley-belt  system  to  a  DC-motor.  The  variable  speed  motor  is 
calibrated  such  that  a  desired  constant  angular  velocity  in  either  CW  or  CCW  direction  can  be 
maintained.  To  insure  rigid-body  rotation  of  the  gas  within  the  tube,  three  cylindrical  layers  of 
honeycomb,  each  1  cm  thick  with  2  mm  hexagonal  holes,  are  successively  positioned  within  the  tube 
near  the  middle.  The  honeycombs  also  help  to  prevent  flame  flash-back  into  the  inlet  supply 
chamber.  Butane  (commercial  grade)  and  air  are  respectively  metered  by  soap-bubble  meter  and 
rotameter  and  fully  mixed  prior  to  introduction  into  the  burner. 

For  lean  flames,  which  have  closed  tips  (Fig.  la)  it  is  found  that  the  rotation  first  results  n  the 
buckling  of  the  central  region  of  the  flame  cone  as  shown  in  Fig.  lb.  This  flame  corresponds  to  the 
fuel  concentration  (percent  by  volume),  the  mean  velocity  over  the  tube  cross  section  =  2.69 
(percent  by  volume),  the  mean  velocity  over  the  tube  cross  section  V  =  122  cm/s,  and  the  burner 
angular  velocity  co  =  5 1  RPS  (rotations  per  second).  At  larger  co,  the  flame  cone  begins  to  enter  the 
tube  (Fig.  lc)  and  eventually  at  co  =  69  RPS  the  flame  flashes  back  into  the  tube  and  assumes  the 
inverted  conical  shape  shown  in  Fig.  Id.  The  inverted  flame  cone  is  convex  towards  the  fresh  gas 
with  peripheries  bending  towards  the  tube  walls,  such  that  the  flame  surface  near  the  walls  is  almost 
concave  towards  the  fresh  gas.  When  similar  experiments  are  performed  with  a  larger  diameter 
tube,  20  mm,  the  transition  from  buckled  state  (Fig.  lb)  to  flash-back  (Fig.  Id)  occurs  very  rapidly.  In 
such  cases,  direct  photography  of  the  intermediate  state  such  as  shown  in  Fig.  lc  will  be  quite 
difficult. 

The  tip  of  the  initial  rich  butane/air  flames  Qf  =  5.60,  V  =  129  cm/s  and  co  =  0,  's  open  [6-8]  (Fig. 
2a).  As  co  is  increased  to  44  RPS,  the  height  of  the  original  flame  cone  is  reduced  and  a  new  inverted 
conical  flame  surface,  convex  towards  the  fresh  gas,  appears  in  the  middle  of  the  outer  curtain-flame 
(Fig.  2b).  The  inner  conical  flame  is  separated  from  the  outer  flame  by  a  dark  ring  of  extinction 
region.  As  co  is  further  increased,  the  inner  conical  flame  begins  to  enter  the  tube  (Fig.  lc). 
Eventually,  at  co  =  50  RPS,  the  inner  conical  flame  flashes  back  into  the  tube  while  the  outer  flame 


cone  vanishes  (Fig.  2d).  Now,  as  opposed  to  the  lean  flame  (Fig.  Id),  the  peripheries  of  the  conical 
flame  are  nearly  parallel  to  the  tube  wall. 


IH.  MATHEMATICAL  MODEL 

In  a  recent  study  of  the  effects  of  centrifugal  and  Coriolis  acceleration  on  the  form  and  stability  of 
a  premixed  flame  freely  propagating  in  a  rotating  tube,  the  following  evolution  equation  for  the 
asymptotic  flame-front  has  been  derived  [51, 


f  +  2  V"  (f  f  =  D'J  2  I  !r  f ) +  4  Y“2  Vi'  (f2 ' 2  ^ 


(3.1) 


where  -z  =  F(r,t)  is  the  flame  front  surface;  vn  is  the  normal  velocity  of  flame  propagation;  is  the 
thermal  diffusivity  of  the  gas;  0)  is  the  angular  velocity  of  rotation;  R  is  the  radius  of  the  tube;  y  is  the 
thermal  expansion  coefficient  (y  =  (pu  -  pb)/Pb);  Pu,  Pb  *  densities  of  the  unbumed  and  burned  gas. 

If  the  tube  is  thermally  insulated,  Eq.  (3  1)  should  be  solved  subject  to  the  boundary  condition 


It  (R’°  =  0 

dr 


(3.2) 


Equation  (3.1)  corresponds  to  the  limit  of  weak  thermal  expansion  (y  «  1)  and  slow  rotation 
(o)R/vn  «  1).  For  simplicity,  the  Lewis  number  is  taken  to  be  unity,  since  the  basic  features  of  the 
flame  response  to  rotation  (buckling  and  flash-lock)  do  not  seem  to  be  very  sensitive  to  the 
composition  of  the  mixture  (Figs.  1  and  2).  Equation  (3.1)  is  written  in  the  frame  of  reference  in 
which  the  velocity  of  the  oncoming  gas  flow  v.  is  equal  to  the  normal  flame  speed  vn,  i.e.,  v.  =  vn. 
For  the  subsequent  analysis  it  is  convenient  to  consider  the  flame  front  dynamics  in  the  frame  of 
reference  where  v.  *  vn. 

The  corresponding  evolution  equation  is  obtained  from  Eq.  (3.1)  via  simple  transformation 


F  =  F  -  (v.o  -  vn)t 


(3.3) 


Hence,  the  modified  evolution  equation  reads 


f  +  2  V"(f  F  =  D,h  r  |(r  f  K  Y“2  Vi‘  (r2  -  2  R2) +  V- '  V" 


(3.4) 


Replacing  the  adiabatic  boundary  condition  (3.2)  by 


F(R,t)  =  0 


(3.5) 


one  obtains  the  system  describing  an  axisymmetric  flame  held  inside  the  rotating  tube,  when  the 
axial  velocity  of  the  oncoming  flow  is  v.  (Fig.  3a).  With  all  its  dissimilarity  to  the  Bunsen  burner 
(Fig.  3b)  this  system  proves  capable  of  capturing  the  main  ingredients  of  the  phenomenon  observed 
in  rotating  Bunsen  flames.  Introducing  new  nondimensional  variables 


(3.6) 


the  problem  (3.4),  (3.5),  may  be  written  as 


at  2 lap/ 


=  1-2- 
P  3p 


a  (  ao 

p 


aP 


+  x(p2 


(3.7) 


<D(l,t)  =0 


(3.8) 


where 


1  _  YC02R4  M  _  (v°»  *  VnjVnR 
A>  —  f  Li.  —  ■ 

4D§,  Dl 


(3.9) 


IV.  TRANSFORMATION  TO  LINEAR  PROBLEM 

For  the  further  analysis  it  is  convenient  to  make  the  following  transformation 
<D=-2  In)  u) 

which  converts  the  nonlinear  problem  (3.7),  (3.8)  to  the  linear  one: 


(4.1) 
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To  exclude  singular  solutions  one  has  to  impose  the  additional  condition: 


~  (0,t)  =  0 

3p 


(4.4) 


The  time-independent  solution  u(°)(p)  of  the  problem  (4.2),  (4.3),  and  (4.4)  corresponds  to  a  steady 
configuration  of  the  Bunsen  flame  d>(o)(p).  For  the  subsequent  arguments  it  is  convenient  to  convert 
the  problem  (4.2-3)  as 


3u 

3s 


(4.5) 


^(0,s)  =  0  ,  u^*^,  s)=l 


(4.6) 


where 

q  =  4fkn  p  ,  s  =  im  x ,  (o  <  %  <  4fm) 

and 


L  =  lf?72-^l2A  (Fig.  4) 
2  2 


Let  us  consider  an  auxilliary  time-independent  problem 


A  (O) 

v<°>(0)  =  1  ,  (0)  =  0  ,  0  <  %  <  oo 


Thus, 


v(°)(5)  =  u(°)(^)/u(°)(0)  (4  9) 

For  sufficiently  small  L  (e.g.  L  <  0)  v(°)(^)  does  not  have  zeros  at  %  >  0  [9].  The  zeros  appear  only 
at  sufficiently  large  values  of  L.  The  first  zero  £i(L)  of  v(°)(^>  emerges  from  +»  at  L  exceeding 
some  critical  value  Li  and  then  monotonically  decreases  to  zero  as  L  — *  «» (Fig.  5).  The  point  ^  = 
corresponds  to  the  singularity  of  the  function 


=  -2  lnju(°l|  .  (4.10) 

For  sufficiently  small  L,  the  function  <l)(o,(^)  has  no  singularities  in  the  interval  0  <  c,  <  4 V  X/2  . 
Since  3L/9X  >  0  (Fig.  4),  the  interval  (0,  4YX/2)  expands  as  L  increases.  Therefore,  when  X  exceeds 
some  critical  value  Xi  the  first  zero  (singularity)  %  =  will  enter  the  interval  (0,  ATkjl).  With 
further  increase  of  the  parameter  X,  the  second  singularity  ^2  >  will  enter  the  interval  (0,  4 V  X/2 ) , 
then  the  third  £3  >  ^2.  and  so  on. 


V.  BASIC  PROPERTIES  OF  THE  EQUILIBRIUM  SOLUTION 

When  p  is  of  order  of  unity,  the  transition  from  one  equilibrium  configuration  <j)(°  to  another 
occurs  quite  gradually.  However,  in  typical  laboratory  situations  4  is  a  rather  large  parameter. 
Indeed,  let  vn  =  25  cm/s,  v*,  =  1.5  vn,  R  =  1  cm,  Dth  =  0.25  cm2/s.  Then,  according  to  (3.9)  p  =  5  x 
103.  In  this  case  the  buckling  (i.e.  change  of  sign  of  <t>pp(0))  occurs  at  X  =  2p  =  104.  In  this  domain 
of  parameters,  even  very  small  variations  of  p  and  X  causes  a  drastic  change  in  the  flame  front 
shape.  This  phenomena  may  be  easily  explained.  Since  d^j/dL  — » <*>  as  L  Li  (Fig.  5),  large  qi  is 
very  sensitive  to  small  variations  in  L  (i.e.  X  and  (t).  On  the  other  hand,  the  singular  solution 

emerges  at  =  4TX/2,  which  is  a  gradually  increasing  function  of  X.  As  a  result,  transition  from  a 
slightly  budded  flameconfiguration  to  a  singular  one  at  sufficiently  large  X  occurs  in  an  almost  jump 
like  manner. 

Such  strong  parameter  dependence  (stiffness)  of  the  rotating  flames  makes  numerical  (and 
laboratory)  study  of  the  system  a  rather  difficult  task.  Figures  (6,7)  show  some  samples  of  the 
numerical  experiments  which  we  have  performed.  Figure  6  corresponds  to  flame  configurations  at  X 
=  32  and  different  values  of  L  (i.e.  p.).  The  buckling  occurs  at  L  =  0  (p.  =  16),  the  first  singularity  at  L 
-  2.25  (p  ~  -2),  i.e.  when  the  flow  in  the  burner  is  lower  than  vn  (3.9). 

Figures  7(a,b,c)  show  flame  shapes  at  X  =  2  x  104.  The  buckling  emerges  at  L  =  0  (p=  104). 
Transitions  from  the  buckled  state  to  the  singular  one  occurs  in  virtually  a  discontinuous  fashion 
somewhere  between  L  =  3.4935286  and  L  =  j. 4935287  (Fig.  7b).  Double  precision  Runge-Kutta 
method  was  use^  for  the  numerical  calculations,  providing  up  to  20  significant  figures,  with  local  error 


of  about  10' 10.  The  integration  was  initiated  from  R  =  1  to  avoid  the  singular  behavior  at  the  axis  of 
symmetry.  In  Figure  7c  (L  =  20,  4=  6  x  103)  flame  front  exhibits  several  singular  points. 

As  L  increase  singularities  shown  in  Figures  7(b,c)  actually  appear  first  at  the  wall.  The 
configuration  represented  by  the  curve  7(b)  corresponds  to  the  situation  in  which  the  singularity  is 
already  shifted  from  the  wail  towards  the  axis  of  the  tube. 


VI.  STABILITY 

It  is  not  difficult  to  show  that  nonsingular  configurations  (X  <  A.i(p.),  Sec.  IV)  are  stable  and 
singular  ones  (A.  <  Ai)  are  not  stable.  Indeed,  let  us  put 

u(£,s)  =  u(°)(£)  +  w(£,s)  (6.1) 

where  disturbance  w(^,s)  is  a  solution  of  the  problem 


dw  _  ^  d 

dT~$d$ 


(6.2) 


w^(0,s)  =  0  ,  w{4VA/2 ,  s)  =  0 


w(^,s)  may  be  sought  in  the  form 


w(5,s)  =  X  wn(^)ea"s 

n=l 

One  obtains  the  following  eigenvalue  problem  for  wn(^) 


;  jL(t  dwn 

%  <$>  r 

1  +  (a„  -  i;2)wn  =  o  , 

where 

o 

ii 

s 

.  c 

.  w„(0)  =  0 

An  —  L  -  on 

,  L  =  -^  fk/2  -  —  flfk 
2  2 

(6.3) 


(6.4) 


(6.5) 
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It  is  known  from  the  general  Sturm-Liouville  Theory  [9],  wn(£,)  has  (n-1)  zeros  in  the  interval 
(0,  4VX72)  and  An+i  >  An  (or  an+i  <  an).  Hence,  as  follows  from  the  results  of  Sec.  IV, 

Ai  =  L  -  01  >  X  or  01  <  0  for  X  <  Xi  (6.6) 

Thus,  an  <  0  for  any  n,  and  therefore  at  X  <  Xi  the  equilibrium  solutions  are  stable.  For  the  same 
reason  at  X  >  Xi 

Ai  =  L  -  0i  <  L  or  0i  >  0  (6.7) 

Therefore,  at  X  >  Xi  the  equilibrium  solution  will  be  unstable  at  least  to  the  disturbances  correspond¬ 
ing  to  the  first  eigenfunction  wi(£).  In  other  words,  all  the  equilibrium  configurations  involving 
singularities  (Fig.  7c)  are  unstable. 

Let  us  now  look  at  the  nonlinear  evolution  of  the  perturbed  solution  d>(°)(£)  in  the  case  of 
instability 

oo 

<D(^,s)  =  -21n|  u(4,s)|  =  -21n|  u<°>(4)  +  £  wn(^)e^  . 

n=l  (6.8) 

Since  On  >  0n-fl» 

0($,s)  ~  *21n  |  u(°)(^)  +  wi(^)e0is  I 

as  s  — >  «>.  Hence,  outside  of  the  boundary  layer  (near  %  =  4VX/2,  (6.9)  yields 
0(i;,s)  ~  -20js  as  s  -»  oo 

Thus,  at  X  >  Xi  as  a  result  of  instability  flame  flashes  back  into  the  burner  tube, 
by  the  burner  rim  is  possible  only  at  sufficiently  slow  rotation. 


VII.  CONCLUDING  REMARKS 

The  proposed  mathematical  model  shows  that  the  observed  buckled  flames  (Figs.  1,2)  emerge 
not  as  a  result  of  an  instability  of  the  classical  Bunsen  cone,  but  as  the  only  possible  equilibrium 
configuration.  The  pertinent  smooth  solutions,  however,  exist  only  at  sufficiently  slow  rotations. 
When  the  angular  velocity  to  exceeds  a  critical  value  ©c  'mooth  solutions  are  replaced  by  oscillating 
singular  solutions  (Fig.  7c).  The  singular  solutions  (unlike  the  smooth  ones)  are,  however, 
unstable.  Thus,  at  ©  >  ©c  and  fixed  velocity  of  the  oncoming  flow  (v.),  flame  stabilization  in  a 
rotating  tube  becomes  impossible  and  the  flame  flashes  back. 


(6.9) 

(6.10) 

Flame  stabilization 


The  physical  reason  for  the  flashback  is  quite  simple.  If  co  >  wc  the  rate  of  fuel  supply  (v.)  turns 
out  to  be  insufficient  to  provide  equilibrium  flame  configuration  in  the  gas  such  that  the  flame  breakes 
away  towards  the  oncoming  flow  of  the  fresh  mixture. 

To  verify  this  interpretation  in  a  more  formal  way  it  is  instructive  to  consider  flame  stabilization 
in  a  wide  tube,  where  the  dissipative  effects  (Dih)  may  be  ignored.  The  corresponding  mathematical 
problem  is  thus  reduced  to  consideration  of  the  first  order  differential  equation 


2VM  =  47“2v"1r2'2RV 


v„  -  Vn 


(7.1) 


Equation  (7.1),  as  is  readily  seen,  admits  a  real-valued  solution  (defined  over  the  entire  cross 
section  of  the  tube,  0  <  r  <  R)  only  at 


to  <  cog  =£V2 v„rl(v„  -  vn) 
K. 


(7.2) 


It  is  interesting  that  at  small  nonzero  Dth,  the  real-valued  equilibrium  solutions  of  Eq.  (3.4)  survive 
at  to  >  tog  due  to  the  'tunnel  effect',  i.e.,  emergence  of  high  frequency  singular  solutions.  These 
solutions  are,  however,  dismissed  via  stability  arguments  (Sec.  VI). 


Acknowledgement 

This  research  was  supported  by  Air  Force  Astronautics  Laboratory  under  Contract  No.  F04611- 
87-K0067  and  DOE  Grant  DE-FG02-88ER13822. 


REFERENCES 

1.  Bidinger,  D.  F.  and  Sohrab,  S.  H.  (1985)  unpublished. 

2.  Chen,  Z.  H.,  Liu,  G.  E.,  and  Sohrab,  S.  H.  (1987)  Combust.  Sci.  and  Tech.  51,  39-50. 

3.  Lin,  T.  H.  and  Sohrab,  S.  H.  (1987)  Combust.  Sci.  and  Tech.  52,  73-79. 

4.  Sivashinsky,  G.  I.  and  Sohrab,  S.  H.  (1987)  Combust.  Sci.  and  Tech.  53,  67-74. 

5.  Sivashinsky,  G.  I.,  Rakib,  Z.,  Matalon,  M.,  and  Sohrab,  S.  H.  (1987)  Combust.  Sci.  and  Tech.,  to 
appear. 

6.  Sivashinsky,  G.  I.  (1975)  J.  Chem.  Phv.  62(2),  638-643. 

7.  Buckmaster.  J.  D.  (19791  Combust.  Sci.  and  Tech.  20.  33. 

8.  Mizomoto,  M.,  Asaka,  Y.,  Ikai,  S.,  and  Law,  C.  K.  (1984)  19th  Symposium  (Int.l  on  Combustion. 
The  Combustion  Institute,  Pittsburgh,  p.  1933. 

9.  Courant,  R.  and  Hilbert,  D.  (1953)  Methods  of  Mathematical  Physics.  John  Wiley,  New  York, 
Vol.  1. 


98 


Figure  Captions 


Figure  i. 

Figure  2. 

Figure  3. 

Figure  4. 
Figure  5. 
Figure  6. 
Figure  7a. 
Figure  7b. 

Figure  7c. 


Direct  photographs  of  lean  butane/air  Bunsen  flames  with  £2p  =  2.69  and  (a)  to  =  0.  (b) 
to  =  51  RPS  (c)  0)  =  69  RPS  (d)  to  =  69  RPS. 

Direct  photographs  of  rich  butane/air  Bunsen  flames  (a)  Qp  =  5.6,  to  =  0  (b)  Qp  =  5.6,  to 
=  44  RPS  (c)  Qf  =  5.6,  to  =  50  RPS  (d)  Q?  =  7.12,  co  =  67  RPS. 

Diagram  illustrating  Bunsen  flames  stabilized  inside  the  tube  (a)  and  at  the  rim  of  the 
burner  (b)  under  absence  of  rotation. 

The  ^.-dependence  of  the  parameter  L. 

The  L-dependence  of  the  First  zero  of  the  function 

Flame  configurations  at  X  =  32  and  different  values  of  the  parameter  L. 

Stable  flame  configurations  with  X  =  2  x  104  and  L  <  Lq  =  3.49357287. 

Transition  from  (1)  stable,  L  =  3.49357286  to  (2)  unstable,  L  =  3.49357287  flame 
configuration  at  X  =  2  x  104. 

Unstable  flame  with  multiple  singular  points  with  X  =  2  \  104,  L  =  20. 
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Figure  4. 
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ABSTRACT 


Premixed  flames  stabilized  in  stagnation  point  flow  are  considered.  It  is  shown 
that  this  classical  flow  field  is  capable  of  sustaining,  apart  from  the  known  planar  flames, 
also  some  nonplanar  flame  configurations  similar  to  those  recently  observed  in 
experiments. 
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1.  INTRODUCTION 


Stagnation  point  flow  is  one  of  the  most  convenient  hydrodynamic  systems  for 
experimental  and  theoretical  studies  of  premixed  laminar  flames.  The  remarkable  feature 
of  the  stagnation  point  flow  is  that  being  essentially  nonuniform  it  is  capable  of  sustaining 
ideally  planar  flame  fronts.  However,  recent  experimental  observations  of  Ishizuka, 
Miyasaka  and  Law  (1982),  Ishizuka  and  Law  (1982)  have  shown  that  the  stagnation 
point  flow  may  stabilize  also  some  nonplanar  flame  configurations.  One  of  these  new 
possibilities  of  the  flame  equilibrium  is  shown  on  Figure  1. 

The  present  paper  is  aimed  to  show  that  the  nonplanar  equilibrium  flame 
configurations  are  indeed  admissible  by  the  classical  stagnation  point  flow  field. 


2.  MATHEMATICAL  MODEL 


Consider  the  idealized  scheme  of  Figure  2  in  which  the  axisymmetric  jet  of  the 
stagnation  point  flow  field  is  assumed  to  be  infinitely  wide  in  diameter  but  of  finite  length 
(h).  The  coordinate  frame  is  chosen  in  such  a  way  that  z  =  0  cooresponds  to  the 
stagnation  plane.  Assuming  for  simplicity  the  gas  to  be  incompressible,  the 
corresponding  axisymmetric  flow  field  v  =  (vr,  vz)  is  given  by 


(2.1) 

v*-vffl[,+(2)] 


where  V  is  the  velocity  of  the  oncoming  gas  flow  (at  z  =  -  h.  Figure  2). 


Define  the  flame  as  a  geometrical  surface,  F(r,z,t)  =  0,  propagating  relative  to 
the  gas  flow  with  a  prescribed  speed  V0. 

The  flame  front  evolution  equation  may  then  be  written  as 


v-  n  -  D  =  -  V„ 


(2.2) 


where 


is  the  flame  velocity  relative  to  the  fixed  frame  of  coordinates  (r,z)  and 


(2.3) 
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(2.4) 


-  VF 
n  =  - 


M 


is  the  normal  to  the  flame  front. 

To  avoid  dealing  with  implicit  functions  let  us  specify  F(r,z,t)  in  two  different 


ways 


and 


F(r,z,t)  =  z  - 
F(r,z,t)  =  r  -  V(r,t; 


In  the  first  case 


-*  _  (dv,  - 1 ) 


n  = 


Vi  +<p? 


,  D  = 


0: 


Vi  +  <£ 


while  in  the  second 


n  = 


,  D  =  - 


Vt 


Vi  +  VTT^ 

Hence  the  evolution  equation  (2.2)  may  be  written  either  as 


vr<fc-vt+  =  -  V0  V  1  +  <fc2 
(z  =0(r,t) 


or  as 


-Vr  +  VWz  +  Vt  =  -  V0V 1  +  \|£ 
(r  =V(z,t) 


Substituting  (2.1),  Eqs.  (2.9)  and  (2.10)  yield 

idi+8ft+^)(2+3+^-v"/7^  ■ 

■  i l )[ 1 1 + si  •  Mh)  [2 + ©1 * + '*'• = • 


(2.5) 

(2.6) 


(2.7) 


(2.8) 


(2.9) 


(2.10) 


(2.11) 

(2.12) 
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Introduce  the  nondimensional  quantities 


p  =  z/h  ,  £  =  z/h  ,  O  =  <yh  ,  4/  =  \|//h  , 

x  =  V0t/h  ,  Q  =  V/V0  >  1 

In  terms  of  the  new  variables  Eqs.  (2.11)  and  (2.12)  become 

<DX  +  V  1  +  <Dp  +  Qp(l  +  cb)  Op  +  Q3>  (2  +  <d)  =  0 
N't  +  V  1  +\|^  -  Qy(l  - 4)  -  C£(2  +  5)  V4  =  0 


(2.13) 


(2.14) 

(2.15) 


3.  EQUILIBRIUM  CONFIGURATIONS  OF  THE  FLAME  FRONT 

From  now  on  we  confine  ourself  to  the  study  of  the  time-independent  flame 
configurations  only,  i.e.,  we  set  =  0.  Equation  (2.14)  then  yields 

^  =  Q2a>2(2  +  <P)2  -  1 _ 

V  Q202(2  +  O)2  +  Q2p2(  1  +  <l>)2  -  1  -  0*0(1  +  <t>)(2  +  <D)p  (3' 1 } 

Here  we  have  selected  the  branch  for  which  d>p  is  negative  at  q>  =  0. 

Figure  3  shows  the  corresponding  integral  curves.  Note  -that  in  the  region 
marked  by  (#)  (characterized  by  relatively  low  velocity  of  gas  flow)  flame  stabilization  is 
impossible.  The  flame  configuration  (Fig.  4b)  corresponding  to  the  solution  (a)  is  of 
special  interest  since,  despite  its  rather  unusual  appearance,  it  seems  to  be  realizable 
experimentally  in  counterflow  systems  (Ishizuka  and  Law  1982).  Similarly,  the  time- 
independent  version  of  Eq.  (2.5)  yields 


Q2(1+4)V-1 _ 

Vq2(1  +S)2  +  Q^2(2+$)-  l  -0^(1  +  5)(2+5)v 


(3.2) 


This  equation  corresponds  to  the  branch  for  which  is  positive  at  t,  =  0.  The  integral 
curves  of  Eq.  (3.2)  are  shown  on  Fig.  5.  In  this  case  each  flame  front  configuration  seems 
to  be  experimentally  feasable,  provided  the  flame  is  held  at  the  rim  of  the  burner 

{%  =  -  1,  \j/  =  \|/0).  Curve  (a)  clearly  corresponds  to  the  Bunsen  type  flame,  while  (b)  to 
the  "grammophone''  flame  (Fig.  4d,c).  Since  Eqs.  (3.1)  and  (3.2)  are  nothing  but  different 
branches  of  the  same  basic  equation  (2.2)  their  solutions  may  well  coexist.  In  other 
words  actual  flame  front  shape  may  be  composed  from  the  elements  generated  by  Eq. 

(3.1)  as  well  as  by  Eq.  (2.2). 


Figure  4e  shows  one  of  such  configurations,  which  clearly  corresponds  to  the 
experimental  situation  shown  on  Fig.  1. 


4.  CONCLUDING  REMARK 

The  simple  geometrical  model  discussed  above,  shows  that  the  classical 
stagnation  point  flow  field  apart  from  the  known  planar  flame  is  indeed  capable  of 
sustaining  other  more  complex  flame  configurations.  In  the  present  article  we  didn't 
discuss  the  question  of  dynamic  stabilization,  which  due  to  apparent  nonuniqueness  of  the 
possible  equilibrium  states  may  prove  to  be  quite  nontrivial. 
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Nozzle  stabilized  counter  flow  lean  propane-air  flame.  (Courtesy  of  S.  H. 
Sohrab,  Northwestern  University.) 

Diagram  illustrating  premixed  flame  moving  in  stagnation-point  flow  Field. 
Shadowed  area  represents  the  burnt  gas  region. 

Integral  curves  of  Eq.  (3.1)  for  <p  =  2.  (a  =  1/Q,  P  =  -  1  +  V 1  -  1/Q  ) 

Integral  curves  of  Eq.  (3.2)  for  cp  =  2.  (a  =  1/Q,  p  =  -  1  +  V 1  -  1/Q  ) 

Some  of  the  possible  equilibrium  flame  configurations  in  counter  (stagnation 
point)  flow  Fields. 
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